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CENTRAL LIMIT THEOREMS FOR THE SPECTRA 
OF CLASSES OF RANDOM FRACTALS 

PHILIPPE H. A. CHARMOY, DAVID A. CROYDON AND BEN M. HAMBLY 


Abstract. We discuss the spectral asymptotics of some open subsets of the real line with 
random fractal boundary and of a random fractal, the continuum random tree. In the case 
of open subsets with random fractal boundary we establish the existence of the second order 
term in the asymptotics almost surely and then determine when there will be a central limit 
theorem which captures the fluctuations around this limit. We will show examples from a 
class of random fractals generated from Dirichlet distributions as this is a relatively simple 
setting in which there are sets where there will and will not be a central limit theorem. 
The Brownian continuum random tree can also be viewed as a random fractal generated by 
a Dirichlet distribution. The first order term in the spectral asymptotics is known almost 
surely and here we show that there is a central limit theorem describing the fluctuations 
about this, though the positivity of the variance arising in the central limit theorem is left 
open. In both cases these fractals can be described through a general Crump-Mode-Jagers 
branching process and we exploit this connection to establish our central limit theorems 
for the higher order terms in the spectral asymptotics. Our main tool is a central limit 
theorem for such general branching processes which we prove under conditions which are 
weaker than those previously known. 

MSC: 28A80 (primary), 60J80, 35P20 (secondary). 


1. Introduction 

Let D be a non-empty bounded open subset of R d for d > 1 and let A be the Dirichlet 
Laplacian on D. Then the spectrum A of —A is discrete and forms a positive increasing 
sequence 

0 < Ai < A 2 < • • • , 

where the eigenvalues are repeated according to their multiplicity. Interest in the geometric 
information about D encoded by A started a little over 100 years ago and was crystallised 
by Kac in his paper |29j entitled ‘Can one hear the shape of a drum?’ Or more precisely, 
does A determine D up to isometry? The answer to that question is no in general, as shown 
in [2D ®j; see also [9j for a concise presentation of a family of counterexamples. 

However some geometric information about D can be recovered. Weyl’s theorem shows 
that the eigenvalue counting function N defined by 

N{\) = #{V : V < A} 

has asymptotic expansion 

N( A) = Cl {d)wo\ d {D)\ d / 2 + o(\ d / 2 ), 

as A -A 00, for some constant c\(d) depending only on d, where vol^ denotes the d-dimensional 
Lebesgue measure. Aside from prompting Kac’s question this result has led to a large body 
of work on the behaviour of the eigenvalue counting function and we now give a very brief 
description of the results that have motivated the work we will present here. 

As a first extension it is natural to ask about the second order term in this expansion. If 
dD is smooth, then under some assumptions, that there are not too many periodic geodesics, 
the expansion has a second order term 

N{ A) = Cl (d)vol d (D)X d / 2 - c 2 (d)vol d -i(dD)\^- 1 ^ 2 + o( A (d “ 1)/2 ), 
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as A -A oo, for some other constant 02 (d) depending only on d. The reader is referred to 
B51 12105103 EDI and references therein for more information. This means that, under some 
regularity conditions, we can recover the size of the domain and that of the boundary from 
the spectral asymptotics; in particular, using the isoperimetric inequality, we can determine 
whether or not D is an open ball. 

Interest in the second term of the expansion of N grew further when Berry studied the 
spectral asymptotics of domains with a fractal boundary in mm- He conjectured that the 
Hausdorff dimension of dD should drive the second order term. Brossard and Carmona in 
[ 8 ] studied the associated partition function, a smoothed version of the eigenvalue counting 
function, and showed that the Minkowski dimension, dM , was the relevant notion of dimen¬ 
sion for the second order term in the short time expansion of this function. For the counting 
function itself a general result of Lapidus [35j shows that, if rf — 1 < < d , the second 

order term is of order 0(\ dM / 2 ) provided the Minkowski content of the boundary is finite. 
In general it is difficult to determine the precise order of growth for the second order term 
for arbitrary boundaries, however for one-dimensional domains [38] it was shown that the 
Minkowski dimension captures the order of growth of the second term in the asymptotics 
and the Minkowski content, the constant, when they exist. 

The problem of determining the spectral asymptotics has also been considered for sets 
which are themselves fractal. For some classes of fractal, such as the Sierpinski gasket, or 
more generally p.c.f. self-similar sets [32] or generalised Sierpinski carpets [3], a Laplacian 
can be defined and shown to have a discrete spectrum. The exponent for the leading order 
growth rate in the eigenvalue counting function is called the spectral dimension and differs 
from the Hausdorff or Minkowski dimension of the set. If the fractal has enough symmetry, 
such as for instance the Sierpinski gasket, then a Weyl type theorem is no longer true 021 , 
[5] in that the rescaled limit of the eigenvalue counting function does not converge. However 
the Weyl limit does exist for ‘generic’ deterministic p.c.f. self-similar sets [33] and also for 
random Sierpinski gaskets [23] and it is natural to ask about the growth of the second order 
term in these settings. 

Our aim is to consider some random fractals where we anticipate more generic behaviour 
of the counting function. We will consider both domains with fractal boundaries and fractal 
sets here. Firstly we will consider the case of open subsets with fractal boundaries in the 
one-dimensional case of a so called fractal string. Our second case will be an example where 
the set itself is a fractal, the continuum random tree. In both cases the first order terms in 
the spectral asymptotics due to the fractal structure are understood and we will focus on 
the behaviour of the second order terms. 

A fractal string is a set obtained as the complement of a Cantor set in the unit interval, so 
can be thought of as a sequence of intervals of decreasing length m- The Dirichlet Laplacian 
is then the union of the Dirichlet Laplacians on each interval. Some discussion of the spectral 
asymptotics of random fractal strings can be found in m where it is shown that for Cantor 
sets constructed via random iterated function systems, the second order term due to the 
boundary exists almost surely. We will consider a suitable subset of these random fractal 
strings and determine when the order of the fluctuations about the boundary term is given 
by a central limit theorem (CLT). 

This turns out to be a subtle question and the existence of a CLT is determined by the 
rate of convergence in an associated renewal theorem. We will give a precise statement after 
introducing all the terminology in Theorem |4.3| We will then show that when the fractal is 


generated using a Dirichlet distribution, the existence of a central limit theorem depends on 
the particular Dirichlet distribution considered. 

An example of what we are able to show is the following. Let S 7jQ; , for 7 G (0,1), a G N, be 
the random fractal string obtained as the complement of the random Cantor set generated 
by subdividing any interval of length £ into three, retaining two intervals of size £, T^ 1 £, 

and removing one of length £(1 — — T^ 7 ), where the pair (Tf, T 2 ) is independent for each 
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interval and distributed as Dirichlet(<a, a) (that is a Beta(<r,<a) distribution in this simple 
case) and 0 < 7 < 1. We write P for the probability law for the random fractal string and E 
for expectation with respect to P. We note that 7 will be the Minkowski dimension of the 
random Cantor set P-almost surely, that is the dimension of the boundary of the string. We 
write N^ a (X) for the associated eigenvalue counting function. 

Theorem 1.1. (i) For all a £ N and 7 G (0,1) there is a strictly positive deterministic 
constant C( 7, a) such that as X -A 00 

A _7//2 A 1 / 2 — W 7 , a (A)) -A (7(7,0) P -almost surely . 

(ii) If a < 59, then there exists a strictly positive deterministic constant a (a) such that as 
A —y 00 

A 7 / 4 ^A -7 / 2 A 1 / 2 — iV 7 jQ! (A)^ — ( 7 ( 7 , 0 )^ -A Z, in distribution 

where Z is normally distributed with mean 0 and variance <j(a ) 2 G (0, 00). 

(Hi) There exists an a > 80 and a 7 G (0,1) such that: if 59 < a < a, then there exists a 
not-identically-zero periodic function p lja (x) such that 

EiV 7 a (A) = 1 a 1 / 2 - C( 7 , a )\^ 2 +p 7 ,a(log A)A 7 ^/ 2 + o(A^), 

7r 

where 17 (a) = max{9ft(0o) G (—00, 1 ) : P(0q) — 0 }, 

pw -nV+o+i)-^ 

2=0 

and, for this range of a we have 1/2 < 7 ( 0 ) < 1 . In particular 

A t/4 (W 2 ( 1 a 1 / 2 - N^ a (\) - C(7,«)) 

does not converge in distribution as A -A 00. 

Remark 1 . 2 . ( 1 ) The first result gives the almost sure behaviour of the second term in 
the counting function asymptotics and is true for random fractal strings constructed using 
a wide class of distributions on the simplex. 

(2) In part (iii) we conjecture that it is possible to take a = 00 and any 7 G (0,1). Indeed, 
towards proving the above result, we first provide conditions under which a CLT holds 
(see Theorem |4.3| and Section [ 5 T 2 ] ) , and explain when one will not (see Remark |4.4[ ). This 
distinction is determined by the rate of convergence in a related renewal theorem and depends 
on the values of the roots of P(9) = 0 , which we solve numerically (we can also solve this 
equation analytically for small values of a). These computations demonstrate that we can 
take a to be at least 81. Furthermore, although we are not able to prove it rigorously, the 
monotonicity of the results suggests that a can be taken arbitrarily large. 

(3) We also conjecture that, in the case where there is no CLT, i.e. a > 59, the size of the 
second order term is determined by 7 ( 0 /), in that, P almost surely for e > 0, 

jv (A) = 1 a 1/2 - C( 7 ,a)A 7/2 + o(A 7??( “ )/2+e ), 

’ 7F 

where 1/2 < rj(a) < 1 and 17 (a) -A 1 as a -A 00. 

(4) The proof of the above result shows that the period of p 7 ^ a is given by 47r/7|3(0o )| 5 where 
60 is one of the complex conjugate pair of roots whose real part gives 17 (a). 

Observe that, as a increases, the Beta(cq a) distribution becomes closer to the distribution 
given by a delta measure at the point (1/2,1/2). If we take 7 = In2/ ln3, then we anticipate 
that our random fractal string should converge (in a suitable sense) to the Cantor string (the 
string formed as the complement of the classical ternary Cantor set) as a goes to infinity. It 
is known that for the Cantor string there is a non-constant periodic function that appears 
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in the second order term in the counting function asymptotics m ■ Thus our result suggests 
that there is a non-trivial transition in the parameter space from the case where there is 
‘enough randomness’ for a CLT about the second order term, to the case where there is not, 
through to the limit, where there is not even a strong law of large numbers for this term. 

We will also consider the case of the Brownian continuum random tree, a random self¬ 
similar fractal. It was shown in mi that there was a Weyl limit for the counting function 
in this case. It was also shown that the second order term for this fractal set was of order 
1 in mean - which would be anticipated as the boundary of the tree is just two points, 
a O-dimensional set. In this paper we show that there is a CLT about the almost sure 
asymptotics. However at this point we have not shown strict positivity of the variance due 
to the complexity of the correlation structure in the variance of the limit of the rescaled 
counting function. We conjecture that there will be a non-trivial CLT for this counting 
function. This will show that the randomness in the construction means the second order 
term in the spectral asymptotics is determined by the fluctuations about the leading order 
term, as these are much greater than the effects due to the boundary of the set. 

The main technical tool we develop is a central limit theorem for the general Crump-Mode- 
Jagers branching process. In our setting the random fractal sets, the random Cantor set 
boundary of the string, or the continuum random tree, can be encoded as general branching 
processes. We are able to use a characteristic associated with these processes to determine 
the behaviour of the counting function. In this case there may be dependence on the offspring 
of an individual and we obtain a CLT in this more general setting, extending the work of [28] . 
We also remark that the techniques used here can easily be applied to geometric counting 
functions or other functions associated with heat flow, such as the partition function or heat 
content of the set. We anticipate similar behaviour in the fluctuations of these quantities 
about their almost sure limits. 

The paper is organised as follows. In Section [2] we recall the definition of the general 
branching process and some laws of large numbers for such processes. We then prove our 
central limit theorem for the general branching process using a Taylor expansion proof. In 
Section [3j we restrict ourselves to general branching processes where a suitable function of 
the birth times is chosen to lie on an n-dimensional simplex, which will ensure that the 
limit of the usual branching process martingale is a constant. We will call such processes 
A n -GBPs and discuss extensively how to establish the conditions required for the central 
limit theorem in this setting as this will allow us to illustrate when we do and do not have a 
central limit theorem for the associated general branching process. In Section [4j we define a 
family of open subsets U of [0,1] whose random boundary is a statistically self-similar Cantor 
set built using scale factors on the simplex. We are then able to show our main result which 
gives conditions for the existence of a central limit theorem. In Section 5 we consider some 
examples where the law of the A n -GBP is given by a Dirichlet distribution. We show that, 
for some Dirichlet weights, the eigenvalue counting function of the set U satisfies a central 
limit theorem. As a consequence we will be able to establish Theorem 0 In Section 6 we 
turn to the continuum random tree. We recall that this tree can be viewed as a random self¬ 
similar set and how to construct a Laplace operator on it. We then show that the conditions 
for the general branching process central limit theorem hold and hence there is a CLT in the 
spectral asymptotics. 

Notation. For convenience, we will use the shorthand notation C{ with i G N to mean some 
positive constant whose value is fixed for the length of a proof or a subsection. 

2. A CENTRAL LIMIT THEOREM FOR GENERAL BRANCHING PROCESSES 

2.1. General branching processes. In this subsection, we introduce the general or C-M-J 
branching process. The presentation is inspired by USESHE to which the reader is referred 
for further information. 
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In the general branching process, the typical individual x is born at time a x , has offspring 
whose birth times are determined by a point process £ x on ( 0 , oo), a lifetime modelled as a 
non-negative random variable L X: and a (possibly random) cadlag function <f> x on R called a 
characteristic . 

We index the individuals of the general branching process using the address space 
( 1 ) / = (J where N° = 0 . 

k> 0 

The ancestor 0 is born at time a 0 = 0, and individual x has £^(0, 00 ) offspring whose birth 
times cr Xj i satisfy 


£ 4 °°) 



2=1 

where S is the Dirac measure and x,i is the concatenation of x and i. The trace of the 
underlying Galton-Watson process is a random subtree of / which we denote by E. We write 
dE for the set of infinite lines of descent in the process. For x, y E E we also use the notation 
x < y if there exists a sequence ( 27 ,..., z\f) with Zi E N, i = 1,..., k with k E N such that 
y — (x, 27 ,..., Zk). Similarly, for x E E, y E dE we write x < y if there exists a sequence 
( 27 , Z 2 , • • •) with 2 ^ E N, i = 1,2,... such that y = (x, 27 , ^ 2 , •..). A cut-set C of E is a 
collection of x E E such that x ^ x' and x' x for all x' 7 ^ x E C and Vy E dE there is an 
t E C such that x < y. 

It is customary to assume that the triples (£ x , L x , 0^)^ are i.i.d. but we allow <f x to depend 
on the progeny of x; we also do not make any assumptions on the joint distribution of 
(£a;, L x: (j) x ). When discussing a generic individual, it is convenient to drop the dependence 
on x and write (£,L,</>). We will write P for the associated probability law and E for its 
expectation. 

We define 

£(£) = £((0, £]), v(dt) = E £(d£), £ 7 (d£) = e _ 7 t £(d£), and z/ 7 (dt) = E£ 7 (d£), 

for 7 E (0, 00 ). Furthermore, we will always assume that the general branching process 
has Malthusian growth , i.e. that there exists a Malthusian parameter 7 E (0, 00 ) for which 
z/ 7 (oo) = 1. This implies, in particular, that the general branching process is super-critical. 
We denote the moments of the probability measure v 1 by 

roo 

( 2 ) Hk= s k u 1 (ds). 

Jo 

In all cases of interest to us, fi 1 will be finite. Note, however, that some convergence results 
still hold when that is not the case, as explained in m ■ 

The presence of the characteristic in the population is captured using the characteristic 
counting process Z^ defined as 

£0(0°) 

( 3 ) z\t) = ^2 <f> x (t - a x ) = <f> 9 (t) + J2 Z t^ ~ 

xCE 2=1 

where the zf are i.i.d. copies of Z^. An important example in the study of random fractals 
is the characteristic </>(£) = (£( 00 ) — £(£))l[ 0jOO )(£), whose corresponding counting process Z^ 
has the property that Z^(t) is the number of offspring born after time t to parents born 
up to time t. Later, we will define characteristics that count eigenvalues of the Dirichlet 
Laplacian. 

There are two central elements in the study of the asymptotics of the counting process. 
The first is that the functions 


z*(t) = e^EZ^t) and u*(t) = e -7 *E <f>(t) 
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satisfy the well-studied renewal equation 

r oo 

(4) Z^{t)=U^(t)+ / Z^(t — S)lSry(ds); 

J 0 

see m for a classic exposition and mmm for alternative results. 

The second is the process defined by 

M t =Y e “ 7 Y 

xeA t 

where 

A t = {x G E : x = (y,i) for some y G E, i G N, and <J y < t < a x } 

is the set of individuals born after time t to parents born up to time t. The process M is a 
non-negative cadlag J^-martingale with unit expectation, where 

Ft = a([F x , a x < t ) and T x = a({(£ y , L y ) : a y < a x }\ : 

we call it the fundamental martingale of the general branching process. 

The martingale convergence theorem shows that M t -A M^ as t -A oo, almost-surely, for 
some random variable M^. Furthermore, under an x\ogx condition standard in the theory 
of branching processes, M is uniformly integrable. More precisely, in cam, Doney proved 
the following result. 

Theorem 2.1 (Doney). The following are equivalent: 

(i) E[£ 7 (oc)(log£ 7 (oo))+] < oo; 

(ii) EMoo > 0 ; 

(hi) EMoo = 1 ; 

(iv) Moo > 0 almost surely on the set where there is no extinction; 

(v) M is uniformly integrable. 

Otherwise, M oo = 0 almost surely. 


For technical reasons, it is often easier to apply renewal theory under the assumption that 
4> vanishes for negative times. When that is not the case, we can set 


£*(oo) 


( 5 ) 

so that 


Xx(t) = c/) x (t)lt>0 + Y Z ti^ ~ £7 *) 1 0 <t«Ji 


2=1 


£ 0 (°°) 


^i[o,oo)W = + Y <Ti> o- 


2=1 


This means that l[ 0 ,oo) — Z x , the counting process of the characteristic y, and we can 
then work with Z x instead of Z^ because x vanishes for negative times and Z x and Z^ 
obviously have the same asymptotics as t -A 00 . 


2 . 2 . Application to statistically self-similar fractals. As discussed in [17, 22,112], the 

general branching process provides a natural way to encode statistically self-similar sets. We 
outline this connection now. 

To build a statistically self-similar set K , we start with the address space / defined in 
0 and a non-empty compact set Kq. To each x € /, we associate a random collection 
{N x , ... $ x ,N x )xei, where N x is a natural number and x ^ are contracting similitudes 
whose ratios we write R x We assume that the collection is i.i.d. in x. 

The random numbers (N x ,x G I) generate a random subtree E of / defined by 0 G E and 

V — Vi ? • • • ^ Un £= S ^—V yi ,..., y n —1 £= ^ and y n Y Ny 1 • 
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For x = x \,..., x n G E, define 

oo 

Kx — ° ° *&x±,...,x n (-^0) nnd iF = |^| 

n=l |x|=n 

where |x| is the length of the word x. The set K has the intuitive property that it can be 
written as scaled i.i.d. copies of itself, namely, 

N 

K=\J<t> i (K i ), 

2=1 

where K \,..., Kjy are i.i.d. copies of K. 

Let us emphasise that the choice of Kq is not unique in general. However, for technical 
reasons discussed in D3, we make the following two assumptions. First, we assume that the 
sets (int K x ,x G E) form a net , i.e. 

x < y =>* int K y C int K x 

and also 

int K x n int K y = 0 if neither x < y nor y < x\ 
the analogue of the open set condition for self-similar sets. Second, we assume that the 
construction of K is proper, i.e. that every cut-set C of E satisfies the condition: for every 
x G C, there exists a point in K x that does not lie in any other K y with y G C. 

The Hausdorff dimension of statistically self-similar sets is almost surely constant on the 
event that it is not empty and was calculated in mmm- It is given in the following result 
by a formula, the random analogue of that due to Moran m and Hutchinson familiar 
from the deterministic setup. 


Theorem 2.2. Let K be a statistically self-similar set. Write (TV, i?i, ..., Rjy) for the num¬ 
ber of similitudes and their ratios. Then, on the event that the set K is not empty, 


dim K — inf < s 


"(I* 1 ) 51 ! 


a.s. 


To specify a general branching process corresponding to the random set K , we set 

N x 

log Rx,n 


2=1 

and L x = sup^cr^ — cr x . With this parametrisation, the set K x in the construction of K 
corresponds to an individual born at time <j x and has size e~ Gx . Furthermore, since 

N \ 

2=1 / 

the Malthusian parameter 7 is equal to the Hausdorff dimension of K by definition. 


E 


e~ sx £(dx) = E 


2.3. Laws of large numbers. Before we can prove our central limit theorem for the general 
branching process, we state Nerman’s laws of large numbers, proved in [45] . They are proved 
for non-negative characteristics with progeny dependence. In applications, if this is not the 
case, it suffices to write the characteristic as the difference of its positive and negative parts. 

We start with the weak law of large numbers. Recall that a measure is said to be lattice 
if its support is contained in a discrete subgroup of R and non-lattice otherwise. 

Theorem 2.3. Let (£ x , L x , <f> x ) x be a general branching process with Malthusian parameter 
7 , where <f> > 0 and <fi(t) = 0 for t < 0 . Assume that u^ is directly Riemann integrable and 
that is non-lattice. Assume further that, for every t, 

sup <f>{u) 

u<t 


E 


< 00. 
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Then, 

roo 

z*(t) -A z^{ oo) = /if 1 / u^(s)ds, 

J 0 

where fi\ is defined in and 

e-^Z*{t) -+ z* (oo)M 00 , in probability , 

as t ^ oo, where M^ is £/ie almost sure limit of the fundamental martingale of the general 
branching process. Furthermore, if M is uniformly integrable, then the convergence also takes 
place in L 1 . 

The strong law of large numbers requires the following additional regularity condition. 

Condition 2.4. There exist non-increasing bounded positive integrable cadlag functions g 
and h on [ 0 , oo) such that 

£ 7 ( qq ) ~ £7(^0 


E 


sup ■ 
.t>0 


9(t) 


< 00 and E 


e ^Mt) 
SUP w.x 

t> 0 h[t) 


< 00 . 


This first part of the condition is satisfied if there exists a non-increasing bounded positive 
function g such that n 7 (l /g) is finite, because then 

£7(00) ~~ £7 ft) 


9(t) 


< 


roo ^ /‘OO 2 

l Wi (i{ds) -L W) (i{ds) ' 


which has finite expectation. As such, this can be thought of as a moment condition that is 
weaker than imposing that z / 7 have a finite second moment; take g(t) = t~ 2 A 1. 

In particular, if the expected number of offspring is finite, this part of the condition is 
satisfied since, with the latter choice of g , 


E 


roo 

/ gft) -1 e ~ l1 £(dt) < sup{(l V t 2 )e _ 7 t }E^(oo) < 00 . 

J 0 t> 0 


We can now state the strong law of large numbers. 

Theorem 2.5. Let (£ x , L x , <f> x ) x be a general branching process with Malthusian parameter 
7 , where <f> > 0 and 0ft) — 0 for t < 0 . Assume that p 7 is non-lattice. Assume further that 
Condition 2Jj_ is satisfied. Then, 

roo 

z*(t) -A 2^(00) = /if 1 / u^(s)ds, 

J 0 

where /i\ is defined at and 

e -7 t Z ( t ) (t) -A z^(oo)M 00 , a.s ., 

as t ^ 00 , where M^ is £/ie almost sure limit of the fundamental martingale of the general 
branching process. Furthermore, if M is uniformly integrable, then the convergence also takes 
place in L 1 . 


Similar results have been proved by Gatzouras in the lattice case. We will not use them 
here and refer the reader to m • 

2.4. The central limit theorem. In [28] . Jagers and Nerman proved a central limit the¬ 
orem for the general branching process under the assumptions that the characteristics are 
i.i.d. We now give a Taylor expansion proof of a similar result, but continue to allow <f x to 
depend on the progeny of x. We start by introducing some additional notation. 

Consider the general branching process (£ X ,L X ,£ X ) X with Malthusian parameter 7 . We 
assume that £ is such that 

Z(t) := Z<{t) 

has zero expectation. In applications, £ is typically a suitably centred version of some 
characteristic 0; we will discuss examples in Section [3] 
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We will use the rescaled version Z oi Z defined by 

f(°°) 

( 6 ) Z(t) = e~^ 2 Z(t) = C 0 (i) + X e-™/ 2 Zi(t - a,), 

2=1 

where £(£) = e _ 7 t / 2 C(£). 

Finally, to have a proxy for the variance, we define 

f(°o) 

(7) y(t) = z(t) 2 = P0 (t) + w* - <*), 

2=1 

where 

£(oo) ?(oo) 

MO = MO 2 + 2 Czi(0 X] - a i) + 2 X I>(*~ 

2=1 2=1 J <2 

As V satisfies an equation of the form ^ which leads to the renewal equation Q, the 
functions v and r defined by 

( 8 ) v(t) = e _7t E V(t) and r(t) — e _7t E p(t) 
satisfy the renewal equation 

POO 

(9) v(t) = r(t) + / v(t — s)v^(ds). 

Jo 

Our central limit theorem requires two technical conditions which we discuss now. 
Condition 2 . 6 . There exists e G (0,1/2) such that 

e -7 */ 2 ( x (t — a x ) -A 0, in probability , 

cr x <et 

as t -A oc. 

This is a regularity condition on (. In applications, we typically expect ( to satisfy a weak 
law of large numbers. Therefore, the sum should grow like e iet and we can expect that the 
condition is satisfied. 

Condition 2.7. There exists n G (0, oo) such that 

supE{|Z(t)| 2+K } < oo. 
t CM 

This is a moment condition. In applications, it is convenient to check it for the third 
moment, i.e. when n — 1 , because that can be done using renewal arguments. 

Theorem 2 . 8 . Let (fi x ,L x , ( x ) x be a general branching process with Malthusian parameter 
7 , where ( is such that E Z(t) = 0 for every t. Assume that v is bounded and that 

v{t) -A v(oo), 


some finite constant, as t -A oo. Assume further that Conditions \276\ and \2.7\ hold. Then, 

m ^ ^oo 5 distribution, 

as t -A oo ; where the distribution of Z ^ is characterised by 


E 




= E 


e -±0 2 v(oo)M Ol 


In the proof, we will use that if z\,... ,z n and w±,... ,w n are complex numbers whose 
modulus is bounded by C, then 


( 10 ) 


ir-n 


2=1 


i ~ 11^2 
2=1 


< c n - x X 


Zi — Wo . 


2=1 


A proof of this may be found in ns Lemma 3.4.3]. 
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Proof. For e G (0,1/2), iterating from <§, the definition of Z, we get 
(11) Z(t)= ^2 e~ 1(Tx/2 Cx(t - a x ) + ^2 e~ lcrx/2 Z x (t - a x ). 

cr x <et xeA et 

The first sum appearing in © can be rewritten as 

e -71/2 ^2 C xit- a x ) 

a x <et 

which, by Condition |2.6[ converges to 0 in probability as t -A oo if we choose e appropriately 
small. For the rest of the proof, we fix such a choice of e. 

We now consider the other sum appearing in ©, and show that it converges in distribu¬ 
tion to Z 0 o as t -A oo. The result will then follow from Slutsky’s lemma. In other words, for 
0 G M, we want to show that 


( 12 ) 


E 


iOY, 


xeA et 


W2 Z(t-a x ) _ g-iflVooJAfoo 


-^ 0 , 


as t -A oo. To do this, write, for an x$ G (0,1) that will be chosen below, 


Ae, t = 


sup | Oe 1 (Tx ^Z x {i — 
xeA et 


Gr 


;)| < X 0 


and split (12) as 

(13) E 

(14) + E 

(15) + E 


■ id 'ExeA et e 7ax/2 Zx(t-a x ) _ e -^ 2 t)(oo)M 00 .^c^ 

J e T,xeA et e-^/ 2 Z x (t-a x ) _ E*eA et e ~ 


;A, t 


o 2 Q 2 v(oo) Ea:SA et e _ g — ^0 2 v(oo)Moo . 


e.t 


We will show that each of these terms converge to 0 as t —>• oo. 
Fix 6 e M and 8 € (0,1). Let xq = xq(S) € (0,1) be such that 


\e z — 1 — z\ < 8\z\ and 


e z - 1 - z - 


< <A| 2 > 


whenever z G C satisfies \z\ < xo. And let r = r(5, 9) G (0, oo) be such that for t > r, 

x~ { 2 +K) \ 0 \ 2 +K supE{\Z(u)\ 2 +K }e-^ Kt ^ < s, 

ueR 

0 2 \\v\\ooe~ iet < x 0 

and 

|p(oo) — p(t/2)| < 5, 


where k is given by Condition 2.7 


Let us start by dealing with (13). For t > r, 


P« t ) = P sup 


xeA et 


|fcr^-/ 2 Z*(( - <t i )| 2+ '‘ > xg+d 


< P ( £ 1*1 

\xeA et 


2 +^ e - 7 <r cc (l+^/ 2 ) 
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which, by Markov’s inequality, is bounded by 


(16) 


x 0 


(2+«) |0|2+K e —7ent/2 


F e 1 <Jx \Z x (t - a x 

.xeA et 


2 T k, 


< x~ { 2 +K) \9\ 2 +K e-^ 2 B 


e-^*E{\Z x {t - a x )\ 2 +K \T £t } 

.xeA et 


< x~ ( 2+k) \9\ 2+k sup E{\Z{u)\ 2+K }e~^ Kt / 2 

uGM 


<5, 


where we have used that a x G T&, that Z x is independent of T e t and that E M et — 1. (This 
is where we need to control the moment of order 2 + k for some k G (0, oo).) Therefore, ( [l3] ) 
is dominated by 

2P OL) < 25- 

Since 5 is arbitrary, it follows that 


E 


e i0J2 xeAet e 


't’ T x/ 2 Z x (t-a x ) _ ~^e 2 v(oo)M 00 . ac 

c i s*e,t 




as t —y 00 , as required. 

To deal with (15), recall that 


e 7<T;r — M e t -A Mqq, a.s., 

xeA et 


as t A oo. Therefore, dominated convergence implies that 


E 


3 A^Voo) e 7<T;z: - e -§ 0 M°°)Mo, 


;A 


e,t 


o, 


as t 


oo, 


as required. 


We will spend the rest of the proof dealing with (14), which we rewrite as 


E 

H E | e ^e-— / 2 Z x (t-a x ) lAx ^ 

V f } 

- p E jeA^ (oo)e 7CT T A?t 

Fet} 


.xeAet 


xeAet 

_ 


using the conditional independence built into the branching structure, where 

Al t = {\ee-^/ 2 Z x (t-a x )\<x 0 }. 


(17) 


By ( |10| ), the term inside the expectation satisfies 

y | E | e iee ~ iax,2 ^x{t~a x ) 


F et }~ I] E {e-^ 2 "( 00 ) e_7CT T^ ( | F tt ] 


xeA et 


1 Af tt |^}-E{ e -^ 2 -(oo)e-w, ^ | j 


xeAet 

< 

xeAet 

A Taylor expansion and the second order exponential estimate yield that 

E <{ e ^e-r^ /2 Z a; (t- ( T a; ) lAa; 




e,t 


-A 


et 


} 


- E 


f 1 + 1 <Tx/2 Z x (t - a x ) - i(9 2 e 1 ,Jx Z x (t - cr,,) 2 ^ 1 a* t 


A, 


et 


~l<7x 
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Furthermore, 


e{(i + iOe 1CTx/2 Z x {t - a x ) - ^(9 2 e 7<T * Z x (t - a x ) 2 ^j 


F 


et 


- E|e"P 2,,(oo)e 7 ax l A : 


< 


E 

+ 


{ (i9e~^/ 2 Z x (t - a x j) l A%t | F et ) 

E { (i - \fe-^Z x (t - a x f - l Alt 


T< 


et 


Now, since 


E 


Z X (t @x) 


F, 


et 


= o, 


reasoning as in (16) shows that 


E { (i0e-^*/ 2 Z x {t - a x )) l A * t | F et ) | 


E { ( i0e-^*/ 2 Z x (t - <r x )) l Q \ <t | F et ) 

< xqEi | XQ 1 \6\e~' 1<Tx ^ 2 \Z x {t - cr x )\l Q \ A x t 

< x 0 E {x~ {2+K) \e\ 2+K e-^ 1+K ^\Z x (t - d x )\ 2+K l n \ Alt \F £t } 

< 

for t > r. And, since for x G A et we have 

\v(oo)6 2 e~ 7(7x /2\ < 9 2 \\v\\ 00 e~ let < xo, 
a Taylor expansion and the first moment exponential estimate yield that 

E { (i - \e 2 e~^Z x {t - a x ) 2 - e —PM°o)e-^ J r £t 


< 


E< 1 ~\° 2e lfTx Z x (t - &x) 2 + ^6 2 v(oo)e 1CTx 


F f 


et 


+ -59 2 \v(oo)\e 1(7x 


< 9 2 {\v(t - a x ) - n(oo)| + 5|M|oo}e J(Tx ■ 
Notice that, for t > r, 

^ 9 2 {\v(t - a x ) - n(oo)| + ^|R||oo}e 


(18) 


xeA et 

<9 2 £ 5(1+11^11^)6-7^+02 £ (2 + 5)11 

x£A et \A t / 2 xeA et nA t/2 

< 69 2 (1 + Halloo) M e t + 9 2 {2 + 5)||p||oo 'y ^ 6 ^ Gx . 

x(zA e tnA t /2 

Together, © to Jl8| ) show that ( [T4l ) is dominated by 


P noC 


~ 1 CF X 


(19) E [5[1 + 0 2 (1 + 2||v|| 00 )]M e t] + E 

xGA e tnA t/2 

Fix c G (0, oo) large. For t large enough, by Lemma 3.5 of [45], 


9 2 {2 + 6 ) !|p||oo 


~7°x 


E 

E 

< E 

E e-7 "- 


x£A e tnA t /2 


£eA et nA et+e 


/ oo 

(1 - v 1 {s))ds 1 
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as t -A oo. The limiting expression can be made arbitrarily small by choosing c large enough. 
Therefore, the second term in (19) converges to 0. Since E M e t — 1, the first term can also 
be made arbitrarily small by adjusting S. The result follows. □ 


2.5. A word on 

that v(t) -A p(oo), 


the lattice case. The statement of Theorem |2.8| requires us to show 
This is typically done using renewal theory. Recall that if z/ 7 is lattice, 


supported on 6Z, say, then this convergence does not occur. Instead we typically get that, 
for every ?/, 

v(y + bn ) -+g{y), 


as n —y oo, where g(y) is a 6-periodic function. In this case, the proof presented above shows 
that, for every ?/, 

Z(y + bn) -A Zoo(y), in distribution, 
as n -A oo, where the distribution of Zoo{y) is characterised by 


E 


e iOZoo (y) 


= E 


e -^0 2 g(y)Moo 


3 . The central limit theorem for A^-general branching processes 

In this section, we discuss applications of the central limit theorem to the general branching 
process under the assumption that the number of offspring £( 00 ) is constant and equal to n, 
say, and that the birth times or,..., cr n are distributed such that, for some fixed constant 
7 E (0, 00 ), we have (e _7cri ,..., e~ 1CTn ) is a distribution on the simplex A n in R n , i.e. 

2=1 

The latter assumption ensures that the fundamental martingale M = 1 . We call a branching 
process with offspring distribution defined on the simplex in this way a A n -general branching 
process. For simplicity, we also suppose throughout this section that v 1 is non-lattice, though 
as discussed in Section |2.5[ we expect subsequential versions of the results below to hold if 
this is not the case. A basic example that we will return to is the case where the distribution 
of (e _7cri ,..., e _7C7n ) is Dirichlet(oq,..., a n ). 

In this case our general branching process L x , 4> x ) x , with Malthusian parameter 7, will 
satisfy a weak law of large numbers, i.e. 

e -7 t z^(t) -a 2 ^( 00 ), in probability, 


as t -A 00 . We wish to describe the random fluctuations around the limit. To do this, we 
study the expression 


( 20 ) 


p 7*/ 2 


e-^Z+it) - z*( 00 ) 


= e“ 7t/2 


Z*(t) - e^z*(t) 1 + e^ 2 [z*(t) - z* (oo)]. 


Our aim is to apply Theorem |2.8| to the first part of this expression. We will see that 
conditions making this possible ensure that the second term converges to 0. This will then 
produce a result on the fluctuations of Z^{t) thanks to Slutsky’s lemma. 

An outline of this section is that we begin by centring our characteristic in order to apply 
the central limit theorem of Section 2. We then show that in order to control the centred 
characteristic we need to control the rate of convergence in the associated renewal theorem. 


3.1. Centring the process. To start, notice that, in the notation of Section [2j 


Z7) - e^z^it) = Z(t) = Zt(t) 
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is a centred characteristic counting process, if we define ( by 

n 

m = Mt ) + E - °i) - ^Z*(t) 


( 21 ) 


2=1 

n 


$$$) + E ( 




2=1 


where we have used that e 7cri H-he 7C7n = 1. This representation can be used to get the 

following control on ( under an assumption on | z^(t — &i )-**(<)!• 

Lemma 3.1. Let (£x,L x ,(/) x ) x be a A n -general branching process. Assume that 

| z^(t) — 2^(00)| < cie _/5lt A C 2 , 
for some positive constants ei,C 2 and (3\ E [ 0 , 7 ]. Then, 

1(001 < 10(01 + 2n ^Cie^ 7- ^ 1 ^ A c 2 e 7 ^ . 

Proof. Notice that 

n 

IC0OI < \Mt)\ + E e 7 (t_<ri) {l^(* - A - A°°)l +1 At) - ^Mll 


2=1 

n 


< \Mt)\ + E (cie (7 " /?l)(t " <Ti) A C2 e 7 ^- <Ti >) + n (cA 7 ”^ 4 A c 2 e 7 

2=1 

/ n n \ 

< |00(OI + ( ci E e (7 “ ft)(t “ CT<) A c 2 E e 7(t_<7<) I + n (cA 7 ^ 4 A c 2 e 74 ) . 


2=1 


2=1 


The result follows upon noticing that 

n n 

e (7-^i)(*-g-») _ e c-Pi)t y~^ g-(7-^i)gj < ne (7-^i)* ; 

2=1 2=1 

and proceeding similarly for the other sum. 


□ 


This lemma and the decomposition in (20) show that understanding the rate of convergence 
of z^(t) to its limit in the renewal theorem is helpful for estimating the fluctuations of ZO 

3.2. Convergence rate in the renewal theorem. Consider the renewal equation 

POO 

( 22 ) z{t) = u(t) + / z{t — s)F{ds ), 

Jo 

where F is a non-lattice distribution function on [0, 00 ). The key to solving the renewal 
equation is the renewal measure 


H = E F* n , 

72=0 

where denotes the n-fold convolution of F with itself, because z is typically given by 

roo 

z(t)= / u(t — y)H(dy); 

Jo 

see da ed] among others. The renewal theorem of m states that if F has a finite mean fi 1 , 
then 

(23) H(t) ~ /xpt, 

where /(x) ^ g(x) means that f{x)/g{x) -A 1 as x -A 00 . Under some regularity conditions 
on ?i, e.g. if u is directly Riemann integrable, this can be used to show that 


1 

:{t) —* z(oo) = — / u(s)ds , 

Mi ./—00 
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as t —» oo; see m- 

The error in the linear approximation of the renewal function 
(24) G(t) = H(t) - upt 


has been studied extensively. When F has a finite second moment /i 2 , 
(25) 




as t — y oo. 

The rate of convergence of G{t) to its limit in this case can be studied using the Fourier 
transform / of F defined by 


/ oo 

e ws F(ds ), w G 

-oo 


We refer the reader to pa sni mi m] and references therein for proofs; see also Appendix B 
of [32] for a discussion of the lattice case. In particular, Stone proved the following theorem 
in 


Theorem 3.2 (Stone). Suppose that there exists r\ G (0, oo) such that f(w) is analytic and 
1 when SRw G (0,rr). Then, for every r G (0,rr) ; 


F2_ 

2 Tl 


G(t ) - ^ = 0(e 


—rt\ 


as t ^ oo. 


Since we aim to apply Lemma |3.1[ we will be particularly interested in exponential rates 
of convergence of z{t) to its limit. This rate of convergence is connected to that in Theorem 


3.2 by the following lemma adapted from H2j which we include for convenience. 

. Suppose that 


Lemma 3.3. Let z, u and F satisfy the renewal equation 

r oo /♦oo 

z{t)= / u(t - y)H(dy) -» up / u(y)dy, 

JO J —oo 

as t — y oo. Then, 

roo roo 

z(oo) - z(t) = tip / u[t + y)dy — / u(t-y)G(dy). 

Jo Jo 

Proof. It suffices to notice that 

roo roc rco 

z(oo) - z(t) = up u(t + y)dy + fip u(t-y)dy- u(t-y)H(dy), 
Jo Jo Jo 


and use the definition of G. 


□ 


3.3. Checking the conditions of the central limit theorem. Discussing the conditions 
of the central limit theorem involves finding growth estimates for £ which can be verified 
using Lemma 3T To be compatible with the framework of Nerman’s law of large numbers, 
we will focus on the situation where the characteristic vanishes for negative times; extensions 


will be discussed where needed. We start by looking at Condition 2.6 


Lemma 3.4. Let (£ X iL x ,(j) x ) x be a A n -general branching process and let ( be defined as in 
(21). Assume that <f>(t) — 0 for t < 0 and that there exists a constant c\ > 0 such that 

ICWI < cie^, 

for some f3\ G (0,7/2). Then, Condition\2.(\ is satisfied. 
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Proof. Let e G ( 0 , 7/2 — /3i). Then, for t > 0, 

; - 7t /2 ^2 at 


- cr i 


a x <et 


< e (ft+e-7/2)t e -et J2 Cl. 

a x <et 


By Theorem 2.3 


—et 


1 -A — / e~ ll dt — - in probability, 

A 7 ^ lJo 


7Mi 


as t -A 0. So the result follows since /?i + e < 7 / 2 . 


□ 


Coupled with Lemma |3.1[ this result hints that we should not expect to have a central 
limit theorem when z^(t) — 2^(00) does not decay at least as fast as e _7t / 2 , the threshold 
for the second term in the right-hand side of (20) to converge to 0. We will explain more 


precisely why this is sharp in Remark 4.4 


Let us now discuss the moment condition for the central limit theorem. The next lemmas 


discuss a set of sufficient assumptions for Condition |2.7| to be satisfied for the A n -GBP. We 
introduce the function 

n 

^(0) = E£y*™ 

2=1 

and note that ^( 1 ) = 1 and that if(Q) is strictly decreasing in 0 . 

Lemma 3.5. Let be a A n -general branching process. Then, we have 


E 

and the sum is finite if y G (1, 00). 
Proof. By monotone convergence, 


£■ 

Lr<EE 


-y i° x 




k =0 


E 

J2 e~ vl ° x 

OO 

J2 e~ yiax 


_.tGE 

k =0 

\x\=k 


Further, conditioning on the birth times, we get that 


E 

J2 e~ yiGx 

= E 

n 

y^ e ~yi° x e -yi{? x ,i-cr x ) 

= *l>(y) E 

J2 e~ y ^ 


\x\=k 


\x\—k- 1 2=1 


\x\=k— 1 


Iterating this and summing over k proves the equality with the infinite sum. Noting that 
ip(y) < 1 for y G (1, 00) completes the proof. □ 

Relying on this, we can produce the following estimate on the third moment of the scaled 
process Z. 

Lemma 3.6. Let (£ x ,L x ,(j) x ) x be a A n -general branching process. Assume that <f>{t) — 0 for 
t < 0, that 

|C(t)| < cie 7t//2 for t > 0 

and that v is bounded. Then, 

|3 


Proof. Notice that 


supE|Z(£)p < 00. 
*>o 


Z{tf = W 9 {t) + '£z i {t-a i ) 3 , 

2=1 


(26) 
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where 

n n 

=(q>(t) 3 + 3(0(£) 2 ^2 Zi(t - <Ji) + 3 (0(t) J2 Zi(t ~ <?i)Zj(t - dj) 


+ ^2 Zi(t — Oi)Zj(t — crj)Zk(t — 07,). 

i,j,k=l 
not all equal 

Therefore, it is clear that \Z{t )\ 3 is bounded by 


(27) 

^2\Cx(t - a x )\ 3 



£<GE 


(28) 

+ 3 ^2 10^ _ 

Tl 

&x)\ y ^ \z X: j{t ~ ^, 2)1 


xGT, 

2=1 

(29) 

+ 3 ^2 Kx(* _ 

Tl 

&x)\ y ^ — ^, 2)1 — ) | 


x£T, 

LJ=1 

(30) 

TL 

+ E E 

^2)1 1 {t <Tj)\ | Z x k(t &k)\ 


sGS i,j,k =1 
not all equal 


To prove the lemma, it is sufficient to check that e~ 3ryt ^ 2 times the expectation of each of 
these terms is bounded, which we do now. 

To deal with (27), note that 

00 

e~W/ 2 EY,\Ut - °x)\ 3 < 4vJ2 e ~ h<7x = 4J2^(3/2) k < 00 , 

xGS zGE k =0 

thanks to Lemma [T5l 

Notice that our assumption on the boundedness of v implies that 

E Z(tf < c 2 e 7 *. 


Therefore, we can control the term corresponding to (28) using that 

e —37*/2 E ^ |^( t _ a x )\ 2 J2 \Zx,i(t - <T Xi i) I 




2=1 


< c?e" 37t/2 E 


e 7 (t-a x ) E [Z X)i (t - a x ,i) 2 \F x } 1 ' 2 




2=1 


< cfc 2 E 


< 


nc x C2 


E^'E* 
2=1 
OO 

E^ 3 / 2 )* 

/c=l 




< OO, 


using that cr^ > cr^ and Lemma 3.5 


To deal with (29), we proceed similarly to get that 


-37t/2 E ^2 |^ x (2 -a x )\ Y2 \ z x,i(t - a *,i)\ I z x,j(t - cr x j ) 
XGS 2,j = l 
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< ne 3 ^/ 2 E ^2 I Cx(t - a x )\ ^ \Z x>i (t - a Xji )\ 2 

XGS 2=1 

oo 

< n 2 cic 2 y] V>(3/2) fc . 

fc=o 


The similar argument needed to bound the term corresponding to (30) relies on the ob¬ 
servation that, if i,j and k are not all equal, then, assuming without loss of generality that 
i is different, we have 


E{| Zi(t - <7i )| | Zj(t - Gj )| |Z fc (2 - (Tk)\\F®} 

= E{| Zi(t - <r<)|| J 0 }E{| Zj(t - Cj )| | Z k (t - a k )\\^} 

< (E [Zi(t - <7 i )V 0 ]) 1 / 2 (E[7(2 - <T i ) 2 |^ 0 ]) 1 / 2 (E[Z fe (2 - <7fc) 2 |A)]) 1/2 - 


Using this and reasoning as above then shows that 

n 

e — 3 7 «/2 E ^ ^ 1^(2 - <T Xi i)\ \Z Xi j(t - a x j)\ I Z x>k (t - <T x>k )\ 

xGS 2 ,j,A:=1 
not all equal 

oo 

< n 3 c 2 y^^(3/2) fe . 

fc=o 


The proof is complete. 


□ 


These results are easily combined to produce a version of the CLT for the case with weights 
on the simplex. 

Theorem 3.7. Let (£ x , L x , 4> x ) x be a A n -general branching process such that <f>(t) = 0 for 
t < 0 . Assume that 

| z^(t) — 2 ^( 00 ) | < cie~^ lt 1 


for some /?i G ( 7 / 2 , 00 ), that 

W)\ < c 2 e /?2t , 

for some /3 2 G ( 0 , 7 / 2 ), and that v is bounded with v(t) -A v(oo) as t -A 00 . Then the CLT 
holds in that 


e ~lt/2 z C ( t ) Zx> , in distribution, 

as t -A 00 , where the distribution of is normal with mean 0 and variance v(oc). 


Proof. It follows from Lemmas |3.1[ 3.4 and |3.6| that Conditions 2J3 and |2.7| are satisfied. 
This, with our other assumptions, gives the required conditions for Theorem 2.8, □ 


4 . Spectrum of random self-similar Cantor strings 

In this section, we discuss the spectral asymptotics of a family of open subsets of [0,1] 
whose boundary is a random self-similar Cantor set generated using a distribution on the 
simplex, called A n -random Cantor strings. Spectral asymptotics for a variety of Cantor 
strings have been studied extensively in [21BSU3S1ES] and references therein. We specialise 
the discussion to A n -random Cantor strings here so that we can study the fluctuations of 
the spectrum using the results of the previous section. 
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Figure 1. First 4 iterations of the construction of K with the distribution 
Dir(l, 1,1) and 7 = 0.6. 


4.1. Construction. Choose 7 E (0,1), n > 2 and consider the random vector (Ti,... , T n ) 
with a probability law on the simplex. Start with the unit interval Kq = [0, 1 ] and replace it 
by n equally spaced intervals of length T-^ 7 ,..., Tn 1 . Replace each of these intervals by n 
intervals created independently with the same procedure. Iterating indefinitely, we obtain a 
decreasing sequence of compact sets {Kj,j > 0), and K = CjKj is a statistically self-similar 
Cantor set. Indeed, in the notation of Subsection | 2 . 2 [ it suffices to set 

R n ) = (n, 77 ..., 7 /7 ); 


the maps (4>i,..., 4>w) can easily be deduced from this. 


The corresponding general branch- 

By 


ing process (<$;,L) (no characteristic just yet) is obtained as described in Subsection 2.2 
construction this is a A n -general branching process and its Malthusian parameter is 7 . Fig¬ 
ure [T| depicts the first 4 iterations in the construction of the set K with the distribution 
Dir(l, 1,1). 

The set in which we are interested in this chapter is U — [0,1] \ K , whose boundary is K 
by construction. Thanks to the Lindelof property, U is a countable union of intervals. As 
such, U is a random string in the sense of [38j ESI EE] and references therein. 


4.2. Fractal dimension. Consider the set U defined above and recall that the Hausdorff 
dimension of dU — K is 7 . Using Theorem |2.5[ we now show that the Minkowski dimension 
of dU exists and is also equal to 7 almost surely. 

Theorem 4.1. The Minkowski dimension of the A n -random Cantor string K is almost 
surely equal to 7 . 

Proof. Consider the characteristic function defined by 

000 = f(00)-£(*)■ 

The corresponding counting process Z^ counts the number of offspring born after time t 
to parents born up to time t. As such, Z^{t) is an upper bound for the covering number 
N(e~ f ^K) of K with balls of radius e ~ l . 

By the strong law of large numbers, we have 

e-^z^it) — / E e-^ s (f)(s)ds € ( 0 ,oo), a.s., 

Ti J 0 

as t —>> 00 . This is easily used to check that diniM^ < 7 almost surely. The result follows 
since dim K — 7 almost surely. □ 


4.3. Spectrum of the Dirichlet Laplacian. Recall that the eigenvalue counting function 
for a domain D (or a countable union of domains) of M is defined by 

Nd( A) = ^{eigenvalues of — A < A}. 


N d ( A) = tvoli(D)A 1/2 - N d ( A). 
7r 


Following m, we define 
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The function Np has the property that if D\ and D 2 are disjoint, then 

(31) N DlUD2 (\) = N Dl (A) + N D2 ( A). 

Furthermore, for r G ( 0 , 00), a change of variables shows that 

(32) N rD ( A) = N D (r 2 A). 

In our applications of the central limit theorem, we will rely on the following assumption, 
which was discussed in the previous section. 


Assumption 4.2. The rate of convergence in the renewal theorem satisfies 


G{t) 


li Z + <T Z 


2 u 2 


< cie 


-pit 


for large t, for some fii G ( 7 / 2 , 00) 

We may now use the general branching process to study Njj. 


Theorem 4.3. Let K be a A n -random Cantor string with dimension 7 and consider the 
string U = [0,1] \ K. Then, 

\- l/ 2 Nu{\) -> 2ft, a.s. and in L 1 , 

as X ^ 00 , for some strictly positive constant 9T. 

Furthermore, if Assumption holds, then 

A 7//4 (A _7//2 7Vt/(A) — Tt) Z in distribution , 

as A ^ oo, where Z has a normal distributions with mean 0 and variance a 2 for some strictly 
positive constant a. 


Proof. Define the random variable S by 

S = ( l~Ri - R n )/(n- 1). 

By construction and the properties in ( piTj ) and (32), we have 

n n 

Nui A) = (n - 1)N S[ 0 , 1 ] (A) + y] AW A) = (n - l)iV [0i i](5 2 A) + ^ N^X), 

2=1 2=1 

where U{ are i.i.d. copies of U. 

Now recall that the eigenvalues of —A for the unit interval [0,1] are (mr) 2 . Therefore, 

%i](A) = 7r _1 A 1/2 - Ltt _1 A 1/2 J, 

which is bounded by 1 A ( 7 t _ 1 A 1 / 2 ). 

To use the general branching process, set 

<£(t) = (n - l)N m (S 2 e 2t ) and Z*(t) = N v (e 2t ) 

so that 

Z^(t) = Ut) + Y J Zf(t-a l ), 

2=1 

where Zf are i.i.d. copies of Z^ and Z^ is the counting process of the characteristic (j). 
Furthermore, notice that 

(33) 0 < <f(t) <{n- l)e t lt < 0 + cil t >o and Z^{t)l t<0 < c 2 e t l t< o, 

for some positive constants c\ and c 2 - 

To establish the first statement of the theorem, we use © and set 


A it) = (p(t)l t >0 + £ i z}{t- CTi) 1 


0 <t<(Ji 


2=1 


which is bounded thanks to (33). 
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Thanks to Theorem |2.5[ this implies that 


r oo 


e lt Z x {t ) —^ 1 / e 7S Ex(s)ds E (0,oo), a.s. and in L 1 , 

Jo 


as t —» oo. By definition, this means that 

roo 

\-^ 2 Nu(\) ^ Ml 1 / ' 

Jo 


>—7 s 


Ex(s)ds, a.s. and in L 1 , 


as A -A oo, as required. 

Let us now prove the second part of the theorem under Assumption T2 Consider and 
£ x defined as in (21) using </> and y. Notice that 

:= z^ x (t) and Z* := (t) 

are equal for t > 0. In particular, the corresponding variance functions v x (t) and v^{t) are 
equal for t > 0. Furthermore, applying the fact that is i.i.d. together with the bounds 

from Lemma 3.1, it follows that there are positive constants c, r = (2 (3\ — 7) A 7 such that 

r*(t) = e-^E\(*(t)\ 2 < ce~ T W. 

These observations and the renewal theorem of m imply that 


(34) 


say. 


roo 

lim v x (t) — lim v^{t) = fiV 1 / e _7t E (^(s) 2 ds = v(oo) £ ( 0 , 00), 
t—> OO t—> OO Jq 


| z x {t) — z x (oo)| < C 2 e ^ 


Since y is bounded and 
(35) 

by Assumption |4.2| and Lemma [3T3| the conditions of Theorem T8 are satisfied. This implies 
that 

e 7t / 2 (e _7t Z x (t) — 91) -A tV(0,p(oo)), in distribution, 

as t 00 . Using the definition of Z x , the decomposition in (20), (35) and Slutsky’s lemma 
completes the proof. □ 

Remark 4.4. The arguments in the proof can be used to show that there cannot be a 
central limit theorem when the rate of convergence in the renewal theorem is not fast enough. 
Suppose that 

z*(t) — 00 ) = cie~ f3lt + o(e _/5lt ), 

as t -A 00 , for some real constant c\. Then notice that the centring introduced in (21) 
and used in the proof of Theorem |4.3| satisfies 

C0 (*) = fa{t) + \z^(t ~ <?i) - z^(oo) + A(00) - z*(t)\ 

2=1 

= <fo(t) + ^ e -(7-/3i)^ _ + o ( e (7-0i)t) 

= ^(t) + aeb-WR + 

say, as t -A 00 (where the remainder is deterministic). Notice that i? is a strictly positive 
random variable. In particular, for some eo £ (0, 00 ), we have P(i? > eo) > 0. Therefore, 
there exists to such that, for t > to, 

IC0 (A 2 > IC0 (i)| 2 lfl>e 0 > 2 C l e o e2(7_/?l)tlii - eo ' 


This implies that, for t > to, 


At) = e-^E|Cj(i)| 2 > \c 2 e 2 e^~ 2 ^{R > e 0 ). 
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If /?i G (0,7/2], then, by renewal theory, we cannot have a finite limiting variance in (34) 
and in particular, no central limit theorem. 


5 . Dirichlet distributions 

We now consider the case where the distribution on the simplex is the Dirichlet distribution 
and perform some explicit calculations to show the range of behaviour that is possible in our 
set up. In this set up the birth times are distributed as 

(e _7<T1 ,..., e~ 1<7n ) ~ Dir(«i,..., a n ), 

where Dir denotes the Dirichlet distribution. We will describe this situation by saying that 
the A n -general branching process has Dirichlet weights a = ( 07 ,..., a n ) and write 

Do = OL\ + • • • + D n , 

as usual. 


5.1. Explicit calculations with Dirichlet weights. In general, it is difficult to determine 
the solutions to f(w ) = 1 needed to use Theorem 3.2 In some cases with Dirichlet weights 
which we discuss now, however, we can use the properties of the Gamma function to study 
/ and deduce convergence rates for the renewal theorem. 

Let X ~ Dir(a) be a random vector in M n . It is well-known that, for every i, 

X{ ~ Beta(c^, Do d^), 

where Beta denotes the Beta distribution. Recall that if Y ~ Beta(/?i,/y then 

-£Y e = + /U/?2) 


where B is the Beta function, i.e. 

Therefore, we get that 

(36) ip(0) = E 


B (x,y) = 


B(/?i,/3 2 ) ’ 

T(x)T(y) 


E x ? 


_ 2 = 1 


r (x + y) ‘ 


T(d 0 ) T(d^ + 6) 


T(d 0 + 9) T(d^) 


where the equation defines -0. For the general branching process with Dirichlet weights 
defined above, it follows that 


/ OO PC 

e ws »~ f (ds) = / 
-OO J — ( 


0 (w-rf)s 


v(ds) = E 


E' 

_ 2=1 


-7(^(1 —W/7) 


= ip(l — w/ 7 ). 


If Do — ol{ G Z for every i, we can use that T(r; + 1 ) = wT(w) to reduce the function if; to 
a rational function which may be simpler to analyse. Notice that this assumption implies in 
particular that there exist some a G M and ^ G Z such that, for every i, we have ol{ — a + ^. 
Therefore, 

Do = cl + £0 — 07 + • • • + OL n — na + t\ + • • • + £ n , 

from which it follows that (n — 1 )a G Z. 

By definition of G in ( |24| ) and writing F = z/ 7 , we have 

r 00 

G * F(t) = H * F(t) - /xp / [t - s)F(ds) 

Jo 

= H(t) — lt>o — 1^1 + 1 

= G{t) + lt<o- 
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Denoting by g the Fourier transform of G and / that of F = z/ 7 , we thus get that 

1 


g(w) = 


1 -f(w)’ 


w G 


Applying (36) and (37), for the general branching process with Dirichlet weights oc satisfying 
(To — G Z for every i, the function / can be written 


n 

/H = y>(i - u>/7) = ^ ~PjTw) 

2=1 ' 


where Pi is a polynomial of degree ao — ol{. It follows that 

nr=i PiH 


g(w) = 




=i+ r” 

= i + 


Er=iii^?> 


nr=i^H-EMRyTiH 

Q(w)' 


say. It is easy to see that 

deg R < (n — l)ao = deg Q. 
Now, decompose g into partial fractions and write 


g{w) = 1 + 


Qi(w) 


^(w- Pi )™i ’ 

where (pi,i < q) are the roots of Q with corresponding multiplicities mi and Qi are polyno¬ 
mials with degQ^ < mi for every i; in particular, 

m\ + • • • + rri q = (n — 1 )<to- 

Recall that, for k G Z + and 3ft (A — r) < 0, 

/ oo p oo 

e x H k e- rt l t > 0 dt = / t k eP~ r ^dt = 

-OO J 0 

and therefore that 


k\ 


(r — X) k+1 


°o d k 




(38) L - ** 

Using this, it is easy to check that 


(i fc e“ rt ) 


dt — 


k\X k 


(r — A) fc+1 ' 


G(dt) = S 0 (t)dt + ^ Qi(t)e Pit l t< odt + y^ Qi,(t)e Pit l t > 0 dt, 

$t Pi < 0 5Rpi>0 


where the Qi are polynomials determined using (38) and satisfying degQ^ < m^. Of course, 
since F is supported on [0, oo), so is H and therefore, by definition of G, we have 


Putting this together shows that 

(39) G(dt) = 5 0 (t)dt - g,^ l l t< odt + y^ Qi(t)e~ pit l t >odt, 

Kpi> o 

which we can integrate to study the asymptotics of G. A particular example which will guide 
us below is given in the following lemma. 
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Lemma 5.1. Assume that all the roots of Q are simple. Then, 

G{t) — — 1 lf<ot + W~2^-t>0 + ~ e P ^l£>Cb 

^>0 pi 

where Ci = Res(g; pf), the residue of g at pi. 

Since in this case none of the singularities can be removed and all have order one, all the q 
are non zero. Furthermore, since if pi is a root with residue q then pi is a root with residue 
Ci as g(w) = g(w), roots with the same real part cannot cancel out. In particular, in this 
case, the result of Theorem |3.2| is sharp. 


Proof. Since all the roots of Q are simple, we have 


g(w) = 1 + 


1=1 


w - pi 


where we must have = Res(g; pi). Integrating ( [39] ) then shows that 

G{t) = lt> 0 - gf l lt<ot - ^ —lt>o+ ^ —e _p ^lt>o- 




Pi 




Pi 


Since F has a second moment, the result now follows from (25) after letting t oo. 


□ 


5.2. Examples. Let us first discuss how the observations above enable us to establish the 
desired rate of convergence for some simple cases of Dirichlet weights. 

Example 1. 

Lemma 5.2. Consider the general branching processes with Dirichlet weights cx described 
above. Assume that 

k 

OR = = a n = - -, k G {1,2, 3 , 4 }, n > 2. 

n — 1 

Then the Fourier transform f(w) of is analytic and ^ 1 when Rw G (0,7]. In particular, 

G(i) - A = O(e-T), 

Z P' 1 

as t -A 00. 


Proof. Letting a — k/(n— 1) a direct calculation gives 

m =n a+ °+!_ i - 9 > 

2=1 

There is always a solution to if (9) = 1 at 9 = 1 and all we require is that the other solutions 
are less than 0 to establish, via ( |37| ), that f(w) is analytic and / 1 on Rw G (0, (1 + 11 ) 7 ). 
For k = 1, the only solution to if (9) = 1 is 9 = 1. 

For k — 2, the other solution to if (9) = 1 is given by 9 — —2 (a + 1). 

For k = 3, the other solutions to if (9) = 1 are 

9 = _lL±i ± V- 3a 2-l2a-8. 

2 2 

and k — 4 has solutions to if ( 6 ) = 1 at 

O 2 

6 = —2d — 4,-d zb — \/ — 4 d 2 — 20d — 15. 

Z z 

Thus the real parts of all the solutions are less than zero and we have the required analyticity. 
The rest of the statement follows from Theorem 13.21 □ 
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Figure 2. Phase plots of 1 — f(jw) for a = 1, 2, 3, 10, 30 and 60. The 
black line indicates the set {z G C : Sftz = 1/2}. The regions of the plot are 
{z G C : Jftz and ^sz G [—8, s]} for 8 = 5, 10, 10, 40, 50 and 50. 


Analytic solutions for the solutions to the equation i/j(9) = 1 do not appear to be available 
for larger values of k. 

Example 2. 


Here we discuss the general branching process derived from the class of examples mentioned 
in the Introduction in which the Dirichlet weights are of the form ol — (cp a) with a G N. 
We will establish the Theorem from the Introduction. 


Thanks to Lemma 5.2, we know that if a < 4, as n = 2, then the Fourier transform / of 


defined in (37) can be used to show that the rate of convergence in the renewal theorem 

In other words, the applicability of 


is sufficiently fast for the requirements of Theorem 2 
Theorem |2.8| depends on the regularity of the characteristic 
More generally we need to solve the equation 


/( t (1 -0)) = ‘ip(9) = 


2r(2a)r(a + 6 ) 


= 1. 


r(a)r(2a + 9) 

As a G N this is a polynomial equation and hence we seek roots of 


a —1 

{0 + OL + l) 

i=0 


2(2<a — 1)! 
O-l)! 


By letting w = 1 — 6 the rate of convergence in the renewal theorem is given by the root 
of 1 — f(jw) with smallest strictly positive real part. We have computed these values 
numerically. 

The numerical evidence shows that when a increases, some roots of 1 — f(jw) get close to 
the imaginary axis. This phenomenon is illustrated in Figure [2] which contains phase plots 
of 1 — /(jw) for different values of cq we rescaled for convenience. To highlight this more 
clearly, Figure [3] contains some close-ups of phase plots showing the absence or presence of 
such roots of 1-/(7^) for different values of a. In particular, when a = 30, the two non-zero 
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Figure 3. Phase plots of l—f(jw) for a = 30, 60, 80. The black line indicates 
the set {z G C : = 1/2}. The region of the plot is {z G C : Jftz G 

[0,2] and G [8,10]}. 


roots of 1 — f(jw) closest to the imaginary axis are 

p± - 0.9951 ±9.1074i; 

when ol — 60, they are 

p± ~ 0.4962 ±9.1027i; 

and when a = 80, they are 

p± ~ 0.3718 + 9.0963L 

We have in fact computed the real part of the relevant root for all values of a from 1 to 80 
- these are plotted in Figure |4j Numerically, this establishes that a = 60 is the smallest 
integer value for which 1 — f(jw) has roots with real part <1/2. 

Our computations also show that for 1 < a < 80 the roots of 1 — f(jw) are all simple and 
occur as complex conjugate pairs except for the root at 0. 

To summarise, this numerical evidence shows that the general branching process with 
Dirichlet weights (cq a) admits a central limit theorem of the type described when a < 59, 
but not when 60 < a < 80. Moreover, the monotonicity of the plot in Figure [4] suggests that 
the range for which there is not a central limit theorem extends to all a > 60. 

We note that we see similar results in the asymmetric case with Dirichlet weights (or, < 22 ), 
or, 0 L 2 C N with ct 2 < cti — 1. In this case the polynomial equation becomes 

/OL2~ 1 / 1 \ , \ Ot-\— Ot2~ 1 


11 (ai + 9 + i) — 


(cki + a 2 — 1)! 
(«i - 1)! 


| [ {pc 2 + 9 + i) — 


(ati + a 2 - 1)! 
(«2 - 1 )! 


\ 2=0 V 7 / 2=0 

Here is a table showing for a given 0 L 2 the values of oq below which we are in the central 
limit theorem regime. 


OL2 

1 

2 

3 

4 

5 

6 

7 

OL\ ~ 1 

OL l 

26 

32 

39 

45 

51 

57 

64 

60 


5.3. Applications to random self-similar strings. For the range of examples considered 

we know that the Cantor set in Figure [l] 


in Example 1 of Section 3, thanks to Lemma 5.2 


satisfies Assumption |4.2| and so, by Theorem |4.3| the corresponding Cantor string satisfies a 
spectral central limit theorem. 

We now return to the second example of Section 3, which was also discussed in the In¬ 
troduction. Figure [5] contains some pictures of statistically self-similar Cantor sets with 
Dirichlet weights (a, a) discussed in Subsection 5.2 The figure illustrates the fact that the 
geometry of the Cantor set becomes more rigid as a increases, because the corresponding 
Dirichlet distribution becomes more concentrated. 


Proof of Theorem M The numerical evidence discussed in Subsection |5.2| shows that As¬ 
sumption [T2] is satisfied for integers a < 59. Thus, by Theorem T3, we have established 
parts (i) and (ii) of the theorem. 
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Figure 4. Plot of the real part of the root of 1 — f(yw) in > 0 closest to 
the imaginary axis against a. 


For (in), we start by noting if S := 1 — T 1 / 7 — (1 — T) 1 / 7 , where T is a [0, l]-valued random 
variable with density 

r(2a )x a - 1 {l-x) a ~ 1 , 


r (a) 

and S := S/n, then the explicit form of yields the following distributional equality: 

(t>{t) = Se l - [Se*}. 


This is clearly bounded above by 1 for all t G 1, and moreover, we recall from (33) that 
4>(t) < e t for t < 0. Taking expectations, the same is true of E 4>{t). Such an observation, 
together with the asymptotic behaviour of the renewal function (as stated at (23)), readily 
allows us to apply the double-sided renewal theorem of m Theorem 5] to deduce that 

roc r oo 

z 4 ’(t)= u^(t — y)H(dy) —>■ /i^ 1 / u (/> (y)dy =: A(oo). 

J 0 J —oo 

Thus we can apply Lemma [T3| to obtain that 

roc i roc 

A(i)-A( oo) = / u <t> (t-y)G(dy) -/ u (t> (t + y)dy. 

Jo Mi Jo 


Using the bounds from ( [33] ) again, it is straightforward to see that the second term is of 
order e _7t . We now examine the first term. Using ( [39] ) , we see 

r oo rt 

/ u*(t-y)G(dy) = / - y)G{dy) 

Jo Jo 

= £ / - y)Qi{y)e-HVdy. 

Kpi> o^° 
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Define (3\ := 7 1 min^ p . > o which by our numerical study in Example 2 of Section 
we know satisfies f3\ G (0,1/2) (for 60 < a < 80). Then 


5.2 


E f 

^Pi>P 1 JQ 


=-7 (t- 


4>(t - y)Qi(y)e PiV dy 


< 


ci Y (1 + t mi ~ 1 )e-^ t [ e^ri y dy 
'a Jo 


®Pi>Pi 




< C 2 J2 (1 + t mi ~ l )e~ lt (l + e 

3lpi>f3i 

= o(e~ Plt ). 

Without loss of generality we label the remaining pair of terms with p± = 7 (/A dzi/^), and we 
have that, as all the roots ar e sim ple and come in conjugate pairs (again, for 60 < a < 80), 
by the remarks after Lemma 5.1, Q±(y ) = ce ±lc for some c, c with c > 0 . Hence 


f e 7 ^ y ^E(j){t — y)Qi(y)e py dy = y^ce ±zc e f e 7 ^ 1 p± ^ v E(j){t — y)dy 
J 0 ± Jo 

= Vce^e-^ f e-^~^ y my)dy. 

1 Jo 


3fy>i=0i 


As 1 —/3i >0 and (j){y ) is a bounded function, the integrals in the above expression converge, 
as t -A 00 , to complex constants Re ±l ° e~ rr ^ 1 ~ p± ^ y E(f)(y)dy. It follows that 

z*(t) — z^{ 00 ) = 2 i?ccos( 7 / 32 t — 0 — c)e~ 1 ^ lt + o(e _7/5lt ). 

Now, if we suppose that R > 0 , then the reasoning in Remark |4.4| indicates that the Cantor 
string does not satisfy a spectral central limit theorem for values of a G {60,..., 80} (recall 
that we have checked numerically that /A < 7/2 and also c > 0 for a in this range). Moreover, 
splitting the process as in (| 20 |) but without scaling, then taking expectations, we can write 


E Z*(t) = e 7 Y(oo) + e 7t (A(t) - ^( 00 )) 

= e yt z^( 00) + 2Rccos(7^2^ — 0 )e 1 ^ 1 ~^ t + o(e 7 ^ 1- ^ 1 ^). 

Rewriting in terms of the counting function we have the result for the mean counting function 
with y(a) = 1 — /3i the required root of the polynomial appearing in the Theorem. 

Thus to complete the proof of (iii) it remains to check that R > 0. We will do this 
numerically for a G {60,..., 80}. First, observe that for a G C with dXa G (0,1), 


/ := 


El e 

'— 00 

00 r ln((n+l)/S) 


= ES" 


= E S a 


/ >00 

e~ at E <f,(t)dt 

-OO 

r e~ at (se f - 77 ) 

J — OO ' 

OO n. 

u 

00 \ 

N a n ) > 


dt 


E > '/ _ e~ at (Se t - n) dt 

Jln(n/S) ' 


(1 + n 1 ) 1 a (1 + n x ) a 
1 — a a 


1 


a(l — a) 


Kn =0 


where ao := (1 — a) 1 and, for n > 1 , 

n 1-a 

x n ” 


07 


a(l — a) 


((1 + n -1 ) -a (l + an -1 ) — l) 
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Figure 5. Statistically self-similar Cantor strings for the distribution 
Dir (a, ck), with a = 1 , 30 and 80 and 7 = 0.6. 


Some elementary complex analysis yields 


(1 + n ) a — 1 + an 


i a{a + 1 ) 


n 


< 16 Mn 3 , Vn > 4, 


where 


M := max 1(1 + z )~ a I < |3a L 

1*1 = 2 


Hence, if n > 4, then 
-l-9ta 


n 


an, 


< 


n 


2 

l-0ta 


— 1 , a ( a + 1) —2\ 


|a(l — a) | 

/|a 2 (a + l) 


(1 — an 1 + n 2 )(l + an 1 ) — 1 — 


a(l — a) 


2 n 2 


+ 16Mn 3 |1 + an 


-l 


|a(l — a) | 

< n _ 2 _ 5 Ra /(a), 

where 


2 n 3 


+ 16Mn 3 11 + an 1 


/(«) := ( |0> „ +1)l + 2‘+»« e SW|l + o|) . 

|a(l a)\ \ 2 J 

Now, 

/»oo /»0 /‘O 

/ e _(rt E</>(f)dt = J - / e~ at E<j){t)dt = I - / e^ES^dt = 1- a 0 ES. 

J 0 J—oo J— oo 

So, setting a = 7(1 — p±), we obtain that for TV > 3, 

/tv N \ 

Re ±id - ES a I ^ a n + §£(1 + a) - § 5T n_(1+a) ~ «o (E5 a - Es) 

\n=l n=l J 


< E^ a 52 

n=A+l 


n 


-l-9ta 


<ES* a 5T n 

n=N+l 


-2-Ota 


/(a) < E5 : 


'!*Ra 


N~ l -^ a f(a) 

1 + fRa ’ 


where £(x) = is f^e usual zeta function. In particular, the above inequality allows 

us to compute an estimate for Re ±l ° whose error is no greater than the upper bound. For 
values of a G {60,..., 80} and 7 = 5 , our computations establish that R > 0, as desired. For 
example, with this choice of 7 , we find that if a = 60, then R ~ 0.09703, and if a = 80, then 
R ~ 0.1056. Note that values of p± and R for all values of a G {60,..., 80} are presented in 
the Appendix below. □ 


6. Spectral central limit theorem for the Brownian CRT 


6.1. Brownian CRT definition and main result. Building on the investigations into 
the spectrum of the Brownian continuum random tree (CRT) undertaken in [11, 12], in this 
section we apply Theorem |2.8| to deduce a central limit theorem for the Brownian CRT’s 


eigenvalue counting function. The starting point for doing this is the characterisation of the 
Brownian CRT as a random self-similar fractal tree with Dir(l/2,1/2,1/2) weights. (This 
was shown in HU using a decomposition first derived in [3]) 
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Figure 6 . An excursion and associated real tree. 

To introduce the Brownian CRT precisely, it will be most convenient to use the now 
well-known connection between real trees and excursions. In particular, a function / is 
said to be an excursion of length l G (0,oo) if it belongs to C(M + ,R + ) and also satisfies 
f(x) > 0 if and only if x G (0,£). Given such a function, define a distance on [0,£] by setting 
df(x,y ) := f(x) + f(y) — 2inf{/(r) : r <E [xAy,x\Zy}}, and let be the equivalence relation 
arrived at by supposing x y if and only if df(x,y) — 0. Subsequently, if 7/ := [0,£|/ 
and dj- f is the corresponding quotient metric, it is possible to check that (7 /,^) is a real 
tree (see m Definition 2.1] for the definition of a real tree, m Theorem 2.1] for a proof of 
this fact, and Figure [ 6 ] for a pictorial example). Applying this construction, one may define 
the Brownian CRT to be the random real tree T — (T, dr) := (72e> d^ e ), where e is simply 
the Brownian excursion normalised to have unit length (see [ 2 j Corollary 22]). 

For P-a.e. realisation of T, it is possible to define naturally an associated measure and 
Dirichlet form as follows. Firstly, the canonical measure on T, which will be denoted by / 17 -, 
is obtained by pushing-forward Lebesgue measure on [0,1] by the quotient map onto 77 This 
procedure yields a non-atomic Borel probability measure of full support, P-a.s. Secondly, 
as a consequence of m Theorem 5.4], it is possible to build a local, regular, conservative 
Dirichlet form (£ 7 -, T'j) on L 2 (T, /py), which is related to the metric dj- through, for every 

dr(x, y)~ l = inf {£ r (f,f) : f £ Tj , f(x) = 0 , f(y ) = 1 }. 

The eigenvalues of the triple (£ 7 -, / 17 -) are defined to be the numbers A which satisfy 

£r(f, 9 ) = A J fgdy T , V# e T r 

for some eigenfunction / G Tj. The corresponding eigenvalue counting function, A 7 -, is 
obtained by setting 

Nj'(X) := ^{eigenvalues of (St^T^t) < 

and it is this function that will be of interest here. We note that it was checked in mi Section 
6 ] that N't is well-defined and finite for any AgI, P-a.s. Moreover, from [Til Theorem 2 ] 
and m Theorem 1.1 and Remark 1.2], we know that there exists a deterministic constant 
Co G ( 0 , 00) such that, as A -A 00, 

(40) FiNj-(X) = C 0 A 2//3 + 0(1), 
and also, P-a.s., 

(41) A“ 2 / 3 lV r (A) ->• C Q . 

These establish second order mean behaviour, and first order almost-sure behaviour of the 
eigenvalue counting function. Here, we further investigate the second order distributional 
behaviour, applying our central limit theorem to prove the following result in particular. 
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Figure 7. Self-similar decomposition of the continuum random tree. 


Theorem 6.1. There exist constants Co G ( 0 ,oo) and C\ G [ 0 ,oo) such that, as X ^ oo, 


in distribution. 


N t (X) - C 0 A 2 / 3 
AV3 




Remark 6.2. Unfortunately we are not able to establish that the asymptotic variance C\ 
is strictly positive, as we were in the corresponding result for fractal strings (Theorem 4-3). 
This is due to the more complicated correlation structure of the relevant characteristics, for 
which we could not find suitable tools to analyse. 


6.2. Self-similarity of the Brownian CRT. As noted above, the key tool in studying 
the spectrum of the Brownian CRT in mm was a self-similar decomposition. We again 
take this recursion as our starting point, and proceed in this section to describe this in more 
detail. We also make the connection with the branching process framework of Section [2| 

Let p G T be the ^-equivalence class of T and x^\ x^ be two ^-random vertices 
of T. Since T is a real tree, there exists a unique branch-point b^(p, x^ l \ x^) G T of 
these three vertices. To be more precise, this is the sole element in the set [[p, x^}} n 
[[x^\x^]] fl [[x( 2 ),p]], where [[x,y]] is the unique injective path from x to y in T. Now, by 
the non-atomicity of /Tp, the vertices p,x^\x^ are distinct almost-surely, and therefore lie 
in different components of T\b r (p, x^\ x^). We will label by 7i, T 2 and 75 the components 
containing p, x W and x^ 2 \ respectively. Moreover, for i = 1, 2, 3, we define a metric dp i and 

_1 /O -1 

probability measure pn on 7} by setting dp := A- dp|p x p, TTi(') : = Ag p(* H7I), where 
A i pj-iTi). Note that, since pp has full-support, A i is almost-surely non-zero. We also fix 
Pi — P 2 — Ps — b T (p,x^\x^), set x^p — p,x^\x^ for i = 1,2,3, respectively, and choose 

xf^ to be a pp-random vertex of 71 for each i — 1,2, 3. (See Figure [ 7 }) A minor adaptation 
of [2J Theorem 2] using the invariance under re-rooting of the Brownian CRT (see [1, Section 
2.7], for example) then yields the following. 

Lemma 6.3. The collections (Ti, djr, pjr, pi, x^\ x^), i = 1,2,3, are independent copies of 
(T, dp, pp, p, x^\ x^), and moreover, the entire family of random variables is independent 
of (A i)f =1 , which has a Dir(^, |) distribution. 


We will label the objects generated by applying this procedure repeatedly using a subset 
of the address space of sequences / introduced in Section [T H In particular, for n > 0, let 
Yi n := {l,2,3} n (using the convention that {1,2,3}° = {0}), and define S := U m >oE m . For 
i G E m , j G E n , we continue to write the convolution ij — i\.. .i m j 1 • •. j n • For k G E, 
we denote by |fc| the unique integer n such that k G E n . We will also write for i G S m , 
i\ n = i\ • • • i n for any n < m. 

Returning to our inductive procedure, given (7i, dp, pp, Pi, where z G E, we 


define (%j, d Tij , UTi,, Pij, %if, ) and A„, j = 1,2,3, from (Ti,d Ti , UTi, Pi,xp ,xf') using 
the same method as that by which T was decomposed above. If the cr-algebra generated by 
the random variables (A^) 1 <|^|< n is denoted by T n for each n G N, then Lemma 6.3 readily 
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yields the following corollary. As in m and m, it is this result that facilitates all that 
follows. 

Corollary 6.4. For eachn G N ; {(71, dp^ pp^ pi, xf\ xf^)}ipz n is an independent collection 
of copies of (F,dp, /ap, p^x^\x^), independent of F n . 


To prove Theorem 6.1, we will work with the Dirichlet eigenvalues of (£77 Fp, fip ). These 
are defined to be the eigenvalues of the triple (£j?, T 7 ^, ^y), where := £\jtd x jtd and 

Fj- := {/ G Fp : f(p) = f(x^) = 0}. Since the corresponding eigenvalue counting function 
(N^( A)) A€ r satisfies 


(42) TVf (A) < Np(X) < iVf (A) + 2 , VA G R, 

(see HD Lemma 19]), the asymptotics of N® are indistinguishable from those of Np at the 
level at which we are working. 

We now make the connection between the eigenvalue counting function on F and a 
general branching process. Suppose that, starting from the single individual 0, each individ¬ 
ual i has three offspring, born at times — |lnA^, j — 1,2,3, after i was born (so that the 
entire population can be indexed by the set £). In particular, this implies that an individual 
i G S has birth time cq — — | lnL^, where D$ := 1 and Di ... A^^ for i G H\{0}. 

For our purposes, we do not need to define lifetimes of individuals explicitly. We do, however, 
define characteristics (<^)^s 5 via the formula 


(43) iv/v) = m +EN/V a3/2 )> 

3 = 1 

where NP is the Dirichlet eigenvalue counting function on (71,^, fip). Note that [TT1 
Lemma 19] implies that <f>i{t) G [0,6] for every t G K, P-a.s. Note also that the random 
function fa only depends on the progeny of i (including the birth times of the offspring of 
i). Thus, we have a general branching process in the sense of Section 2.1, and, in the sense 
of Section [5] it has Dirichlet weights. It is easy to check that this process has Malthusian 
parameter equal to 7 = 2 / 3 . Moreover, iterating ( 43 ) (and checking that the remainder term 
converges to 0 ) allows one to deduce that the corresponding characteristic counting process 

( 44 ) Zp) = J2Mt~P 

satisfies Z^{t) = (see the proof of m Lemma 3 . 5 ]). As before, the rescaled means 

of Z^ and <f will be written z^(t) e _ 7 t E(Z^(t)), u^(t) := e _ 7 ^E(</>(£)), where we omit the 
index from (j> in the expectation since this is unimportant. Both of the above functions are 
well-defined and finite for all f G K (see mi). In fact, 


(45) 


M := supz‘ 

t eM 


V) 


is a finite constant (see [11, Lemma 20]). Moreover, it was proved as [LL, Proposition 21] that 
z*(t) -+ z* (<x>) := u^(t)dt G (0, 00). (The proof that 2^(00) G (0, 00) was actually not 
included there, but this is a simple consequence of flOj Proposition 1.7] and [ID Corollary 
4].) We also have that, P-a.s., e“ 7 t Z^(£) -A z^( 00), see im Proposition 22] - as in the fractal 
strings with Dirichlet weights example, the fundamental martingale is identically equal to 
one, and so the limit is deterministic. Note that a simple reparameterisation of the two 


previous results yields the first order parts of (40) and (41) 


To prove Theorem |6.1[ we introduce a rescaled centred version of the characteristic count¬ 
ing process. Specifically, as before, we set 

Z(t) := Z^{t) = Z*(t) - e^z^it), Z(t) := e~^ t/2 Z{t), 

where £ is defined as at @>. Just as pjj ) was fundamental to demonstrating the first order 
asymptotic behaviour of Np(t) in the arguments of [LL], the recursions at ©> and 0 are 
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central to our efforts to derive the corresponding second order behaviour via the branching 
process result of Theorem 2.8 We note that the use of an analogous recursion formula for 
providing second order bounds was already noticed in [12] • However, that paper was mainly 
focused on the infinite variance D-stable tree case, and did not obtain the type of detailed 
results that we do here for the Brownian CRT. 

6.3. Variance convergence. In this section, we use the renewal equation of Q to show 
that the rescaled variance v{t) := e _7t E (Z(t) 2 ) = E(Z(t) 2 ) converges as t -A oo to a finite 
constant. To do this, we are required to check that v, r and v 1 are suitably well-behaved, 
where r is defined at ^ and i/ 7 (dt) := Xa=i e _7 *P(cq G dt ) - this is the content of the next 
three lemmas. In the proof of the following result, we recall the function fj(x) = 3/(1 + 2x) 
for x > —1/2, as introduced in 


Lemma 6.5. The function v is bounded and measurable, and v(t) —0 as t ^ — oo. 

Proof. We start by checking that v is bounded for t > 0. Similarly to the proof of m Lemma 
5.3], by appealing to [12| Lemma 5.2], it is possible to deduce that v(t) < 2e 7t (/i + I 2 + 73 ), 
where 


h = 


h = 


Y, E ~ Vi) ~ E (Mt ~ *i)|A)) 2 ) , 

E E [ ■ D i\ Y - °ij) - E(A - (nj)\D, 




<3=1 







3=1 


Since <j>(t) G [0,6], I\ can be bounded as follows: 


(46) 


Vies / 


= 6e" 27t E (Z*(t)) = 6 e-^z^it) < 6Me' 7 ', 


where the first equality is a consequence of (|44j), and M is defined as at (|45j). 

For I 2 , first observe that 

3 3 

(47) s 22/\ i j z4 ’(t-a i j) = ^/\ i jz (l> (t-a i j), 

3=1 3=1 

where z^{t) z^(t) — A(oo), and the equality holds because Y^=i Af ~ 1* Now, by results 

of p2, Section 3], we have that \z^{t)\ < Ce~ l1 for t G M. Thus 


A 


< CA 7 * 


^ A ijz^ft - aij ) - E(A^-2^(£ - cTzj ) | D 
3=1 

for some deterministic constant C. In particular, we have proved that 

3 


h < Ce~^J2 E ( Di 


Y ) - E(A ijZ^it - aij)\Di 

3 =1 

Our next step is to show that the above sum is bounded. Writing z^(s,t) := z^(s) — A(t), 
we can proceed similarly to (47) to deduce that 


E 


Y - Vij) - E (DijZ^(t - <Jij)\Di 

3 =1 
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< 2E J ^2 DijZ^(t - aij,t - ci 
\j=i 

From m Section 3], we have for any s < t that = u^(s) — u^{t) — f* u^(w)dw, and 

hence 


E e 

z(EX 


E DijZ^it - cTij) - E (DijZ^it - Oij)\Di 


3 = 1 


(48) 

(49) 

(50) 


< 2j> E D ijU *(t - Gi) 

*se \j =l 


+ 2 ^E 

ies \j =l 


IJ) 


+2 


/ 3 ft-CTi 

E e E D ‘i{ 

zes \j =l Jt-dij 


u^(w)dw 


To bound these expressions, we will apply the following characterisation of z^{t): 
(51) z*(t) = e _7 *E (Z*(tj) = e~^ ]T E (&(t - a)) = ]T E (D iU *(t - <n)). 


zGS 


zGS 


Specifically, the term at (|48j) satisfies 
3 


2^e A DijU^(t -<Ti)\ =2J2 E (A u*(t - <Ti)) = 2z*(t) < 2M. 


zGE \ j=l 




Similarly, the term at (49) is also bounded above by 2 M. Furthermore, the term at (50) can 
be rewritten as 


2 E e (E d A f u ^ t+ 

ytt A - 1 In A' 


w — &i)dw , 


where (A'-)^ =1 is a copy of (Aj)^ =1 , independent of all the other random variables of the 
discussion. Applying (51), this can be evaluated as 


2E ( V A', f 2?{t + w)dw ] < 3ME ( V A' 

\k )~ \U 3 


j ! In A'-1 I < oo. 


Putting these pieces together, we obtain that 

(52) I 2 < Ce~^ 

for some finite constant C. 

Finally, note that 1% satisfies 

h < e -27t ^2 E - a ij)) > 

zee jeE 


(cf. the proof of [1_2, Lemma 5.3]). Again applying (44), the boundedness of 0 and Lemma 
3.51 it follows that 


h < 6e- 2 ^J2 E {DiZt(t-ai 

ZGE 

= 6e“ 7t ^E (Dfz^it-a., 


zGE 
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< 6Me-^^Txp{2) h 


k =0 


(53) = Ce~^\ 

where again M := sup teR (t), and C := 6M/(1 — ^( 2 )) is a finite constant. 

Summing (46), (52) and ([53]), we obtain that v is bounded for t > 0. We now check 


that v(t) is bounded for t < 0 and converges to 0 as t -A — oo. For this, we use the bound 
E {Z^{t) 2 ) < (7e^ 27+e ^ for t G M (cf. [12, Lemma 4.4]), which implies 

v(t) < e ~ 7< (e (V(f) 2 ) + e 27 V(f) 2 ) < Ce 7t (e et + l) . 

Clearly this yields the desired properties of v(t). Finally, to confirm that v is measurable is 
elementary using the fact that Z^(t) is monotone cadlag, P-a.s. □ 

Lemma 6.6. The function r, as defined at &■ is in L 1 (M) and r(t ) —> 0 as |i| —» oo. 

Proof. It follows from the definition of r that, similarly to the proof of Lemma [ 6 ~ 5 | we have 
|r(f)| < 2 e lt (Ji + J 2 + J 3 ), where 

Jr = e" 27 t Var(<^(t)), 

3 


J 2 = Var 


Js = 


,-2-yt 


J2 - a J 

J = 1 


and we will proceed by showing that the statements of the lemma hold for e jt J^, i = 1, 2, 3. 
As in the previous proof, checking the measurability of the functions is elementary, and so 
we will restrict ourselves to finding suitable bounds for them. Firstly, we have 

e ll Ji < e“ 7 *E (0(t) 2 ) < 6 e“ 7t E (0(t)) = 6 u*{t). 

That u^ G L 1 (R) and u^{t) -A 0 as |£| -A oo was established in DU Lemma 20], and so 
the corresponding result for e l1 J\ also holds. For e 7 t J 2 , we consider the cases t < 0 and 
£ > 0 separately. In particular, we have e 7t J 2 < e 7 t M 2 , which clearly demonstrates that 
e 7t J 2 G L 1 ((— 00 , 0]) and e 7t J 2 -A 0 as t -A — 00 . Furthermore, defining z^(t) := z^(t)—z^{ 00 ) 
as in the previous result and recalling once again that \z^{t)\ < (7e _7t , we are able to deduce 
that 

e'Kh = e 7t Var ^ < e 7 * (3Ce~^) 2 = Ce ~A 

which confirms that e 7 t J 2 G L 1 ([0, 00 )) and e ryt J 2 -A 0 as t -A 00 . Finally, for e lf J^ we 
proceed as follows: 


e 7 t J 3 < 3 1 / 2 e“ 7t E(0(f) 2 )E 


< Ce “ 7t/2 ^E(</>(f))E 

< CA(f) 1 / 2 , 



where for the final inequality we use the fact that v is bounded (Lemma 6.5). Now, from the 
proof of HD Lemma 20], it can be seen that (r ^) 1 / 2 G L 1 (M) (and we have already noted 
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that vP(t) -A 0 as \t\ -A oo). Consequently, we have the desired result for e 7t J3. The lemma 
follows. □ 


Lemma 6.7. The measure z/ 7 is a non-atomic Borel probability measure on [0,oo) and also 
/o°° tUj(dt) = 1 . 

Proof. The proof of this lemma is straightforward and omitted. □ 


In view of the preceding three lemmas, the following result is an immediate application of 
the double-sided renewal theorem of m Theorem 5]. 


Proposition 6.8. The function v converges as t -A 00 to the finite constant v(oo) := 
1 -0 o r (t)dt. 

6.4. Verification of Conditions |2.6| and 2.7| . It now only remains for us to check Con¬ 
ditions 243 and |2.7| before we can apply Theorem |2.8| to deduce the desired central limit 
theorem for the eigenvalue counting function of the Brownian CRT. We start by working 

and, to 


towards an estimate for the third moment of Z, which will confirm Condition 
this end, we use another recursion argument. This is similar to the proof of Lemma 3.6 


2.7 


but more involved due to the lack of a uniform bound for Specifically, iterating (26), we 
deduce that for any k E N 

z(tf = Y Wi(t -Vi)+Y 

The following lemma establishes that the expectation of the remainder term here converges 
to 0 as k —y 00 . 


Lemma 6.9. For each t E M, 


lim E 

k—> 00 





= 0. 


Proof. By Cauchy-Schwarz and Lemma |6.5[ 


(54) 



< E (\Z(t)\ (Z*(tf + e 27 *A(f) 2 )) 

< Ce 7t/2 ((E(ZA ) 4 )) 1/2 + e 2lt M 2 


Applying the characterisation of Z^{t) at (44), we have that 


E (Z 4> (t) 4 )= Y E (Mt- Vi)(f>j(t- <Tj)(f> k (t- cr k )(f>i(t- at)). 


Since 

(55) Mt) < 61{ t >-in5i} < 6 e 9lt 5p, 

where Sf is defined to be the diameter of the metric space ( 71 , djf), which is a random variable 
with a finite positive moments of all orders (see proof of [IT], Lemma 20]), it follows that, for 
any 0 , e > 0 , 

E(zA) 4 ) < Ce A6lt Y E ( D i D j D k D i^ S j ls t f Y) 

(56) < Ce Adlt Y E (z?y +e A^ ( 1 +e Af ( 1 +e) L)f ( 1 +€) y /(1+£) . 


Now, suppose E is viewed as a graph tree with edges between i|u_i and i for each i G 
E\{0}, and the subtree of E spanning z, j, fc, l (and the root 0) has shape as shown in Figure 
[ 8 j where we assume that a, 61 , 62 , i, j, A:, l are distinct. It is then straightforward to check 
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Figure 8. A possible configuration of i, j, fc, l. 


from the independence structure of ( Di)i e s that E (D i 
above by 


0 (i+ e ) n^( 1 + e )^^( 1+e ) n^( 1+e ) 




D 


is bounded 


E ( DT (1+e) 


E 


D 


29(1+6) 


n 29(l+e) 

„ 6 l | |a| + l, 


E 


D 


,29(1+6) 
6 2 _ 

n 29(l + e) 

v ^'W + i, 


xE 


jy9( 1+e) 


D 


9(1+6) 


E 


D 


9(1+6) 


Ihl+i, 


D 


,9(l+e) 
Jl |bl 1 + 1 , 


E 


D 


9(1+6) 


D 


9(1+6) 

fc| |9 2 l + l. 


E 


jy 

D' 


9(1+6) 
1 


,9(1+6) 

!| | 6 2 l + l . 


l+b'l+|fc|+KI- 2 IM- 2 IM-4 


which is equal to 

(5?) ^(40(1+ e)) ^ ^(20(1+ e)) ^ l 6 d+l&2|-2|a|-2 ^(0(1+ e )) ^ |i|+li 

where we again recall ^(x) = 3/(1 + 2x ) for x > — 1 / 2 . Since ^(0) < 1 for any 9 > 1 and 
3 e /( 1 +e )'0(6 l (l + e)) 1 /( 1+e ^ -A ^( 0 ) as e -A 0 , if we are given any 9 > 1 , then it is possible 
to choose e > 0 such that 3('0(0(1 + e))/3) 1 / 1+e < 1. By summing (57) over all suitable 
a, &i, & 2 5 h h k-> l f° r such a choice of 9 and e, it follows that the terms ofthe form considered 
contribute at most the finite amount 

/ 1 \ I 3 \ [ 3 \ 


^l-3^^(4<9(l + e))TO 


3*^(2<9(1 + e))^ J \1- 3^^(9(1 + e)) 


to the sum at (56). For other configurations of i,j, fc,/, it is possible to proceed similarly, 


and consequently prove that, for any e > 0, E(Z^(£) 4 ) < Ce^ 47+e ^. 

Returning to (54), the bound of the previous paragraph implies E(|Z(£)| 3 ) < C'e 5 ^/ 2 ( 1 V 
e ct ), and so 


E | | Zi(t - <Tj )| 3 < Ce 5 ^ 2 { 1 V e et )E ^ D 

J \ieYik 

which converges to 0 as k -A oo. 


5/2 


< C^(5/2) k 


□ 


The first main result of this section is the following, which establishes that Condition |2.7 
holds in the present setting. 

Proposition 6.10. We have that sup tGM E(|Z(t)| 3 ) < oo. 

Proof. As a result of the previous lemma, we have that Z(t) 3 — —cq). Hence, from 

the definition of W, we deduce that E(|Z(£)| 3 ) < E (Ki) + E(A 2 ) + E(i^ 3 ) + E(i^ 4 ), where 
K i, K 2 , -A 3 , A 4 are defined to be the terms appearing in equations J27| ) to ( [30] ) respectively, 
and it will be our goal to show that e -t E (Ki) is bounded for i = 1 , 2, 3, 4. 

Applying the bound for <f at (55) and the estimate |£^(£)| = |z^(£) — 2 ^( 00 ) | < Ce~ lf (as 
well as recalling that is a bounded function), it is straightforward to deduce the existence 
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of a deterministic constant C such that, P-a.s., \Ci(t)\ < C( 1 A (e 7t (l + S/))). This bound 
implies \Ci(t)\ — ICi(^)| 1 ^ 2 |Ci(0l 1 ^ 2 — C'e 7t / 2 (1 + 5^ 2 ), and so e~ t 'E(Ki ) is bounded above by 


k =o 


Ce -*E ^ e t ~ rTi (1 + Si) ) = C £ E (a/ 2 ) E(1 + *) = C £ V(3/2) fc , 
V z G S / z G S 

which is finite, because -0(3/2) < 1. 

Secondly, we proceed similarly to obtain that 

AE(K 2 ) < Ce _t E [ £ e^\l + 5]) £ | Zy(t - ay) 


AGS 


3 = 1 


< Ce _t / 3 E ( ^ A E E ((1 + + ^2 + Witt ~ °ij) 

,iGE j=l 


< Ce-^E ( E ^ E E (AA - 

AGS j=l 


1/2 


< 


Ce _t/3 E E(e (t_ <^)/ 3 ) 

AGS j=l 


< C^V’(3/2) fc , 


fc =0 


where the third inequality is a conditional Cauchy-Schwarz estimate (we also apply the fact 
that the moments of Si are finite), and to deduce the fourth we use Lemma [6^5| 

For the third term, we start by observing that, similarly to (54), E(|Z(£)| 3 ) is bounded 
above by 


(|Z(t )| 7/4 (V (/) 5/4 + e 57 t /V(V /4 )) < Ce 77 */ 8 ((E(^(t ) 10 )) V8 + e 57 t / 4 M 5 / 4 ) . 


E 


By making the obvious extensions to the argument applied in the proof of Lemma |6.9[ it 
is possible to check that, for any e > 0, E(Z^(£) 10 ) < (7e( 107+e ^, and hence E(|E(£)| 3 ) < 
Cei 7t /i 2 (e ct Vl). For any a G [0,1], we also have that \Q(t)\ = |Cz(^)HCz(^)| 1_a < C'e( 1_a ) 7t (l+ 
. Putting these bounds together yields 


fi[ 1 ~ a h 


VE (Ks) < Ce~ l E I e (i-a)7(t-«Ti)(i + sl 1 -^) \Zy(t - ay)Z ik (t - a ik )\ 


AG S 


j,k= 1 


< Ce -( 4 -(i-“h)* E ( E (|4-(t - <ry)| 3 |%n) 

jies j,fc=i 


1/3 


xE (| Z ik (t - <T ifc )| 3 |j]j| +1 ) 


1/3 


< Ce“ ( 1 “ ( i“ a) 7 )t (e 7£t V1)E °J 

ViGE / 

b—n \ / 


where the second inequality is an application of Holder (and we bound the Si term similarly 
to how this was controlled when estimating K 2 above). If a — ^, then for t < 0 we obtain 
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from this that e t K% < Cj^kL o ^(3/2)^ < oo. If a = 1, then it is possible to choose e small 
enough so that the above bound implies, for t > 0 , e~ l K^ < C YlkLn ^(17/12) fe < oo. 

Finally, we can proceed as in the proof of Lemma 3.5 to deduce that e ^EijK ^) < 

again. 

□ 


C J2kLo ^(3/2) fc . The additional input needed to do this is provided by Lemma 
This completes the proof of the proposition. 


6.5 


From the proof of the previous result, we have that \Ci(t)\ < C for some deterministic 
constant C. Hence we can deduce Condition |2.6| by applying the same argument as that 


used to establish Lemma [T4| We simply state the conclusion. 
Proposition 6.11. For every e G (0,1/2) ; 

g - 7 t/2 ^2 Q(t - CTj) —> 0, 


( Ti<et 


in probability as t —> oo. 


To complete the proof of Theorem | 6 . 1 [ note that, by dehnition and (42), 

N r (X) - EN t {\) 


AV3 


- Z(lnA) 


< 2A-V3. 


6.10 and 6.11 allow us to apply Theorem 2.8 to deduce the result 


Hence Propositions | 6 l 
with 

/ oo /»oo 

u^{t)dt G ( 0 ,oo), Ci := p(oo) = / z(t)dt G [ 0 , oo) 

-oo J —oo 


Appendix 


The following table contains the approximate valu es of p± and R for different values of a 
with 7 = ^, as required in the proof of Theorem 


1.1 


a 

P± 

R 

59 

0.495347 ± 9.10306i 

0.0964835 

60 

0.503788 ± 9.1027* 

0.0970307 

61 

0.511952 ±9.10235* 

0.0975642 

62 

0.519852 ±9.10199* 

0.0980839 

63 

0.527501 ±9.10164* 

0.0985906 

64 

0.534909 ± 9.1013* 

0.0990848 

65 

0.54209 ± 9.10096* 

0.0995668 

66 

0.549052 ±9.10062* 

0.100037 

67 

0.555805 ±9.10028* 

0.100496 

68 

0.56236 ± 9.09995* 

0.100945 

69 

0.568724 ± 9.09963* 

0.101382 

70 

0.574906 ± 9.09931* 

0.10181 

71 

0.580913 ± 9.09899* 

0.102228 

72 

0.586753 ± 9.09867* 

0.102636 

73 

0.592432 ± 9.09836* 

0.103034 

74 

0.597958 ± 9.09806* 

0.103425 

75 

0.603335 ± 9.09776* 

0.103806 

76 

0.608571± 9.09746* 

0.10418 

77 

0.613671 ± 9.09717* 

0.104545 

78 

0.618639 ± 9.09688* 

0.104902 

79 

0.623482 ± 9.09659* 

0.105252 

80 

0.628203 ± 9.09631* 

0.105594 
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